THE NON-RELATIVISTIC LIMIT OF THE EULER NORDSTROM SYSTEM 
WITH COSMOLOGICAL CONSTANT 
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Abstract. In this paper the author studies the singular limit c — > oo of the family of Euler- 
Nordstrom systems indexed by the parameters k'^ and c (EN^), where > is the cosmological 
- - - constant and c is the speed of light. Using Christodoulou's techniques to generate energy currents, 

QP ^ the author develops Sobolev estimates that show that for initial data belonging to an appropriate 

^rrj , Sobolev space, as c tends to infinity, the solutions to the EN^ system converge uniformly on 

' a spacetime slab [0, T] X M.^ to the solution of the Euler-Poisson system with the cosmological 

' constant k^. 

-t— > ' 

o 

o 

(yQ . 1. Introduction 

CN ■ 

The Euler-Nordstrom system describes the evolution of a relativistic perfect fluid with self- 
interaction mediated by Nordstrom's theory of gravity. In |18j . we introduced the system in di- 

■ mensionless units and showed that the Cauchy problem is locally well-posed in the Sobolev space0 
for > 3. In this article, we study the Newtonian (also known as the "non-relativistic" ) limit 

• of the family of Euler-Nordstrom systems indexed by the parameters k and c (EN^), where is 
I the cosmological constani0 and c is the speed of light. The limit c — > oo is singular because the EN^ 

system is hyperbolic for all finite c, while the limiting system, namely the Euler-Poisson system with 
a cosmological constant (EP^), is not hyperbolic. Using Christodoulou's techniques 6 to generate 
energy currents, together with elementary harmonic analysis, we develop Sobolev estimates and use 
them to study the singular limit c oo. Although we explicitly discuss only the ENJ: system in 
^ ^ this article, the techniques we apply can be generalized under suitable hypotheses to study singular 

I limits of hyperbolic systems that derive from a Lagrangian and that feature a small parameteiQ. 

CO ' Readers who are interested in similar examples of the analysis of singular limits in partial differen- 

04 I tial equations may consult e.g. [4j, [llj . or ^5J. Our main theorem, which we state here loosely, is 

• the following; we state and prove it rigorously as Theorem [H 

OO 

■ Theorem (Non-relativistic Limit). For initial data belonging to the affine Sobolev space 
, with N > A, the corresponding unique solutions to the EN'^ system (with > 0) 

converge uniformly on a spacetime slab [0, T] x to the unique solutions of the EP^^ system 
as the speed of light c tends to infinity. 



(N 
> 



As discussed in [18], we consider the EN^ system to be a mathematical scalar caricature of the 
Euler-Einstein system with cosmological constant (EE^). We now provide some justification for 
this point of view. First of all, like the EE^ system, the EN^ system is a metric theory of gravity 
featuring gravitational waves that propagate along null cones. Second, the main theorem stated 
above shows that if > and c is large, then the EN^ system well-approximates the EP^ system. 
Furthermore, in [15], Oliynyk shows the existence of a class of non-stationary solutions to the EEg 
system that converge to solutions of the EPq system in the Newtonian limit. Hence, for > 0, the 
Newtonian limit of the EN^ system is the EP^ system, while the Newtonian limit of EEq is the EPq 
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^More precisely, we showed local well-posedness in a suitable affine shift of for N > 3, where by "affine shift" 
of we mean the collection of all functions F such that ||_F — V||^iv < oo, where V is a fixed constant array; see 
Section [2] for further discussion of this function space. 

•^The parameter k'^ > is fixed throughout this article. Remark 14.11 contains an explanation of why our proof 
breaks down in the case = 0. 

■^The small parameter is in the case of the ENJ system. 
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system. Based on these considerations, we therefore expec10 that achieving an understanding of the 
evolution of the EN^ system will provide insight into the behavior of the vastly more complicated 
EE^ system. 

1.1. Outline of the Structure of the Paper. 

Before proceeding, we outline the structure of this article. In Section ^ we introduce some 
notation that we use throughout our discussion. In Section [31 we derive the EN^ equations with 
the parameter c and then rewrite the equations using Newtonian state-space variables, a change of 
variables that is essential for comparing the rclativistic system EN^ to the classical system EP^. In 
Section [4l we provide for convenience the ENJ: and EP^ systems in the form used for the remainder 
of the article. From this form, it is clear that formally, limc_>oo EN^ = EP^. In Section [51 we 
introduce standard PDE matrix notation and discuss the Equations of Variation (EOVJ^), which are 
the linearization of the EN^ and EP„ systems. In Section [51 we provide an extension of the Sobolev- 
Moser calculus that is useful for bookkeeping of the powers of c. We also introduce some hypotheses 
on the c— dependence of the equation of state that are sufficient to prove our main theorem. We 
then apply the calculus to the EN^ system by proving several preliminary lemmas that are useful 
in the technical estimates that appear later. Roughly speaking, the lemmas describe the c — > oo 
asymptotics of the ENJ: system. 

In Section [71 we introduce the energy currents that are used to control the Sobolev norms of the 
solutions. One of the essential features of the currents that we use is that they have a positivity 
property that is uniform for all large c. In Section [8l we describe a class of initial data for which 
our main theorem holds, and in Section [9l we smooth the initial data for technical reasons. In 
Section \T0\ we recall the local existence result of |18| and prove an important precursor to our main 
theorem. Namely, we prove that solutions to the EN^ system exist on a common interval of time 
[0, T] for all large c. This proof is separated into two parts. The first part is a continuous induction 
argument based on some technical lemmas. The second part is the proof of these technical lemmas, 
which are a series of energy estimates derived with the aid of the calculus developed in Section [HI 
The two basic tools we use for generating the energy estimates are energy currents and the estimate 
< C ■ ||A<I> — K^$||/^2, for $ G H^. In Section[TTl we state and prove our main theorem. 

2. Remarks on the Notation 

We introduce here some notation that is used throughout this article, some of which is non- 
standard. We assume that the reader is familiar with standard notation for the spaces and the 
Sobolev spaces H''. Unless otherwise stated, the symbols and H'^ refer to LP(]R^) and H'^{M.'^) 
respectively. 

2.1. Notation Regarding Differential Operators. If F is a scalar or finite-dimensional array- 
valued function on R^+^, then DF denotes the array consisting of all first-order spacetime partial 
derivatives (including the partial derivative with respect to time) of every component of F, while 
y(a)^ denotes the array of consisting of all a*'' order spatial partial derivatives of every component 
of F; this should not be confused with V, which represents covariant differentiation. 

2.2. Index Conventions. We adopt Einstein's convention that diagonally repeated Latin indices 
are summed from 1 to 3, while diagonally repeated Greek indices are summed from to 3. Indices 
are raised an lowered using a spacetime metric, which varies according to context. 

2.3. Notation Regarding Norms and Function Spaces. If V is a constant array, we use the 
notation 

(2.3.1) \\FhliE)'=^'\\F-V\\LHE), 



^We temper this expectation by recalling that our proof does not work in the case = and that in contrast to 
the initial value problem studied here, Oliynyk considers the case = with compactly supported data under an 
adiabatic equation of state. This special class of equations of state allows one to make a "Makino" change of variable 
which regularizes the equations and overcomes the singularities that typically occur in the equations in regions where 
the proper energy density vanishes. Furthermore, this change of variables enables one to write the relativistic Euler 
equations in symmetric hyperbolic form. See | 12| and |16l for additional examples of this change of variables in the 
context of various fluid models. 
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and we denote the set of all Lebesgue measurable functions F such that (_e) < oo by L^{E). 

We also define the H^{E) norm of F by 

(2.3.2) \\F\\h^(e)"^{ E \\dsnli(E)Y' ■ 

\a\<N 

Unless we indicate otherwise, we assume that E — 'S? when the set E is not explicitly written. 
If is a map from [0, T] into the normed function space X, we use the notation 

(2.3.3) llli^lll;,,^ sup \\F{t)\\x. 

tG[0,T] 

We also use the notation C'^([0, T],X) to denote the set of fc-times continuously differentiable maps 
from (0,T) into X that, together with their derivatives up to order fc, extend continuously to [0,T]. 

If E' C K*^ {d frequently equals 3,4, or 10 in this article), then C^{E) denotes the set fc— times 
continuously differentiable functions (either scalar or array-valued, depending on context) on E 
with bounded derivatives up to order k that extend continuously to the closure of E. The norm of 
a function F e C^{E) is defined by 

(2.3.4) \F\k,E'^' J2 sup|a5F(z)l, 

where ds represents differentiation with respect to the arguments z oi F (which may be spacetime 
variables or state-space variables, depending on the context). 

2.4. Notation for c— independent Inequalities. If Ac is a quantity that depends on the pa- 
rameter c, and X is a quantity such that Ac < X holds for all large c, then we indicate this by 
writing 

(2.4.1) Ac<X. 

2.5. Notation Regarding Constants. We use the symbol C to denote a generic constant in 
the estimates below which is free to vary from line to line. If the constant depends on quantities 
such as real numbers N, subsets E of M'', functions ^ of the state-space variables, etc., that are 
peripheral to the argument at hand, we sometimes indicate this dependence by writing C{N,E,^), 
etc. We explicitly show the dependence on such quantities when it is (in our judgment) illuminating, 
but we often omit the dependence on such quantities when it overburdens the notation without 
being illuminating. Occasionally, we shall use additional symbols such as Ai,Z, L2, etc., to denote 
constants that play a distinguished role in the discussion. 

3. The Origin of the EN^ System 

In this section, we insert both the speed of light c and Newton's universal gravitational constant 
G into the Euler-Nordstrom system with a cosmological constant and perform a Newtonian change 
of variables, which brings the system into the form (|4.1.ip - (|4.1.8p . A similar analysis for the 
Vlasov- Nordstrom systenH is carried out in [J] . 

3.1. Deriving the Equations with c as a Parameter. 

We assume that spacetime is a four-dimensional Lorentzian manifold M and that furthermore, 
there is a global rectangular (inertial) coordinate system on A4. We use the notation 

(3.1.1) X = {x°,x\x'^,x^) 

to denote the components of a spacetime point x in this fixed coordinate system, and for this 
preferred time-space splitting, we identify t = x^ with time and s = (x^,x'^,x^) with space. Note 
that we are breaking with the usual convention, which is a;" — ct. The components of the Minkowski 
metric and its inverse in the inertial coordinate system are given by 5^,^ = diag(— c^, 1, 1, 1) and 
gi^'^ = diag(— c~^, 1, 1, 1) respectively. We adopt Nordstrom's postulate, namely that the spacetime 
metric g is related to the Minkowski metric by a conformal scaling factor: 

(3.1.2) g^, = e^^g^,. 

^The Vlasov-Nordstrom (VN) model describes a particle density function / on physical space X momentum space 
that evolves due to self-interaction mediated by Nordstrom's theory of gravity. Various aspects of this system are 
studied, for example, in [2], and [3]. 



In (j3.1.2p . (f) is the dimensionlesfH cosmological Nordstrom potential, a scalar quantity. 

We now briefly introduce the notion of a relativistic perfect fluid. Readers may consult [1] or [5] 
for more background. For a perfect fluid model, the components of the energy-momentum-stress 
density tensor (which is commonly called the "energy-momentum tensor" in the literature) of matter 
read 



(3.1.3) 



where p is the proper energy density of the fluid, p is the pressure (this "proper" quantity is defined 
in a local rest frame), and u is the four-velocity, which is subject to the normalization constraint 

(3.1.4) S^.u'^u" = e^^g^^u^u-" = -c\ 
The Euler equations for a perfect fluid are (see e.g. ^) 

(3.1.5) V^T^^^O (z. = 0,l,2,3) 

(3.1.6) V^(nu^) = 0, 

where n is the proper number density and V denotes the covariant derivative induced by the space- 
time metric g. 

Nordstrom's theorjfl fT4] provides the following evolution equatiorQ for cj) : we define an auxiliary 
energy-momentum-stress density tensor 



(3.1.7) 



TIZ = e^'I'Tt"" = c-^e^'^ip + p)u^'u'' + e'^'l'pg^'\ 



and postulate that is a solution to 

(3.1.8) - K^(j} = -47rc-'*Ge4"^trgT = -A-Kc-^Gg^^T^^^ = Airc-'^Ge'^'t' {p - 3p). 
Note that 

(3.1.9) 00 = cfd^.d^^ct) = -c~^dj(t) + Acl) 

is the wave operator on flat spacetime applied to 0. The virtue of the postulate equation (|3.1.8p , as 
we shall see, is that it provides us with continuity equations p.2.7p for an energy-momcntum-stress 
density tensor Q in Minkowski space. 

We introduce the thermodynamic variable rj, the proper entropy density, and close the system by 
supplying an equation of state, which may depend on c. A "physical" equation of state for a perfect 
fluid state satisfies the following criteria (see e.g. [5]): 

1) p > is a function of n > and rj >0. 

2) p > is defined by 



(3.1.10) 



dp 

P^n — 
on 



(3.1.11) 



where the notation |. indicates partial differentiation with • held constant. 
3) A perfect fluid satisfies 

dp 
On 



>0 ^ 
on 



orj 



> with " = " iff 77 = 0. 



As a consequence, we have that cr, the speed of sound in the fluid, is always real for 77 > 



(3.1.12) 



^2 dg ^ dp_ 
dp 



dp/dn\,j 

4) We also demand that the speed of sound is positive and less than the speed of light 
whenever n > and 77 > 0: 

(3.1.13) n>0and?7>0 0<a <c. 



®In l|3.2.12|l . we rescale the dimensionless cosmological Nordstrom potential by multiplying it by c^, and we refer 
to the rescaled potential as the cosmological Nordstrom potential. 

^Norstrom's theory of gravity, although shown to be physically wrong through experiment, was the first metric 
theory of gravitation. 

^Nordstrom considered only the case ft = 0. 
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Postulates 1-3 express the laws of thermodynamics and fundamental thermodynamic assump- 
tions, while postulate 4 ensures that vectors that are causal with respect to the sound cone are 
necessarily causal with respect to the light cone. 

By p.l.lip . we can solve for cr^ and c~^p as c— indexed functions 6^ and 5Hc respectively of rj 
and p : 

(3.1.14) a'^^'&'Mp) 

(3.1.15) c-^p = D\civ,P)- 

Remark 3.1. Wc wiU make use of the following identity imphed by p.l.l2p . (|3.1.14p . and (|3.1.15p : 

(3.1.16) ^{V,P) =&c\v,p)- 

Remark 3.2. Note that c~^p has the dimensions of mass density. As we will see in Section [SJ 
limc_>oo ^ciViP) '^ill be identified with the Newtonian mass density. 

Remark 3.3. We note that the assumptions p > 0,p > together imply that the energy- momentum- 
stress density tensor (j3.1.3p satisfies both the weak energy condition [T^^X^X'^ > holds whenever 
X is timelike and future-directed) and the strong energy condition {[T^^ -l/2g°''^Ta/3gf,u]X^'X'' > 
holds whenever X is timelike and future-directed). Furthermore, if we assume that the equation 
of state is such that p = when p = 0, then (|3.1.12p and (|3.1.13[) guarantee that p < p. It is 
then easy to check that < p < p implies the dominant energy condition (— T'^X'^ is causal and 
future-directed whenever X is future-directed and causal). 

We summarize by stating that the equations fXL^ - (jXTTB)) . ([XTTS)) . p.l.lOp . and (|3.1.15p 

constitute the EN"^ system. 

3.2. A Reformulation of the EN^ System in Newtonian Variables. 

In this section, we reformulate the ENJ: system as a fixed background theory in flat Minkowski 
space and introduce a Newtonian change of state-space variables. The resulting system ()4.1.1|) - 
(I4.1.8P is an equivalent formulation of the EN^ system. We remark that for the remainder of this 
article, all indices are raised and lowered with the Minkowski metric g, so that d^(f> = g^^d^(j). To 
begin, we use the form of the metric p.l.2p to compute that in our inertial coordinate system, the 
continuity equation (|3.1.5p for the energy-momentum-stress density tensor p.l.3p is given by : 

(3.2.1) - a^T^"^ + er^'^a^^ - e-^^V/j^fi^^ {v = o, i, 2, 3), 

where Tj^^^ is given by p.l.7p . For this calculation we made use of the explicit form of the Christoffel 
symbols in our rectangular coordinate system: 

(3.2.2) r^, = 5:d^<j, + 5^9,0 - g^^.g^^dp^'- 
Using the postulated equation (|3.1.8p for 0, (|3.2.ip can be rewritten as 

(3.2.3) = e^^V.T^'' = 8^ [rrx + ^(9^09^ - ]^g_^- 8'^ cj^d^cj, - \g_^^K'<f>^ 

Let us denote the terms from p.2.3p that are inside the square brackets as 0'^". Since the coordinate- 
divergence of Q vanishes, we are provided with local conservation laws in Minkowski space, and we 
regard Q as an energy-momentum-stress density tensor. 

We also introduce the following state-space variables that play a mathematical roleQ in the sequel: 

(3.2.4) i?e = c-V* = e4*$He(?7,p) 

(3.2.5) P^pe^^. 

After we make this change of variables, the components of 8 read 

4 1 1 

(3.2.6) e^-- [Rc + c-^P] e'^u^u" + Pg^" + £^ {d^'(t>d''4> - -g^''d''4>d^^ - -g^'n'^c^''') , 



"physical" quantities are SHc and p. 
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and we replace (|3.1.5p with the equivalent equation 

(3.2.7) a^e^'^^o (1^ = 0,1,2,3). 

We also expand the covariant differentiation from (|3.1.6p in terms of coordinate derivatives and 
the Christoffel symbols p.2.2|l . arriving at the equation 

(3.2.8) 9^(ne'*'^u^) 0. 

Our goal is to obtain the system EN^. in the form (|4.1.ip - (|4.1.8p below. To this end, we project 
p.2.7p onto the orthogonal complemenil^"l of u and in the direction of u. We therefore introduce the 
rank 3 tensor H, which has the following components in our inertial coordinate system: 

(3.2.9) m"' c-'^e^'I'u^'u'' + g^" . 
n is the projection onto the orthogonal complement of u : 

(3.2.10) n'"'uVM = o (1^ = 0,1,2,3). 

We now introduce the following Newtonian change of state-space variable^ 

(3.2.11) v'''^u^/u'' (j = 1,2,3) 

(3.2.12) $ = c20, 

where v = {v^^v'^^v^) is the Newtonian velocity and <f> is the cosmological-Nordstrom potential. 
Relation (|3.2.11|) can be inverted to give 

(3.2.13) ^" = 6-^7^ 

(3.2.14) u^=e-'^7cw^ 
where 

(3.2.15) 7c(v) 



(c2 - |v|2)l/2- 



Remark 3.4. We provide here a brief elaboration on the Newtonian change of variables. Equation 
p. 2. lip provides the standard relationship between the Newtonian velocity v and the four-velocity 
u: if x'^{t) (v = 0, 1, 2, 3) are the rectangular components of a timclike curve in A4 parameterized by 



X 



_ 



t, and T denotes the proper time parameter, then we have that — dtx^ — [dr / dt)-u^ — /u^ 
(j = 1,2,3.) 

Dimensional analysis suggests the approximate identification (for large c) of the cosmological- 
Nordstrom potential <& from (|3.2.12p with the cosmological-Newtonian potential ^cos-Newt^ which 
by definition solves the non-relativistic equation (|4.2.4p below: the cosmological-Newtonian potential 
has the same dimensions as c^, which suggests that when considering the limit c ^ 00, we should 
re-scale the dimensionless cosmological-Nordstrom potential 0, as we did in p.2.12p . Indeed, our 
main result, which is Theorem [4l shows that with an appropriate formulation of the initial value 
problems for the EN^ and EP^; systems, we have that limc^oo ^ — ^cos-Newt- Dimensional analysis 
also suggests the formal identification of i?oo from (|4.2.ip - (|4.2.4p with limc^oo 9^c(??,p) (for now 
assuming that this limit exists), where dldViP) is defined in (|3.1.15p . 

Furthermore, these changes of variables can be justified through a formal expansion (f> — 0(q) -|- 
c~^<^(i) + • • • , Roo = R{o) +c^^R{i) + - • • , in powers of c"^ in equations (|4.1.ip - ()4.1.4p : equating the 
coefficients of powers of c^^ on each side of the equations implies the formal identificationf[3 4>(o) — 
and (A — k2)(/)(;^) = 47rGi?(o)- If we also take into account the equation {A — K'^)^cos-Newt = 'inGR 



■00 



satisfied by the cosmological-Newtonian potential, it follows that $ =^ c^ip w 0(1) « ^cos-Newt- A 
similar analysis for the Vlasov-Nordstrom system is carried out in [4] . 



"'^''We are referring here to the orthogonal complement defined by the Minkowski metric g. 

^^As suggested by Remark 13.21 even though Rc is not a state-space variable, equation 113. 2. 41 1 also represents a 
Newtonian change of variables. 

"'^■^Upon expansion, the formal equation satisfied by 0(q) is (A — k^)(^(o) = 0, and by imposing vanishing boundary 
conditions at infinity, we conclude that (/)(-qj = 0. 
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Upon making the substitutions (|3.2.1ip - (j3.2.12p and lowering an index with the components 
of n in our incrtial coordinate system read (for 1 < j, k < 3): 



(3.2.16) nr"-c-^,^|vp 

(3.2.17) u'^-c-^jy 

(3.2.18) n^o = -7cV 

(3.2.19) Ui'^'c-^^yvk+Si. 
Furthermore, we will also make use of the relation 

(3.2.20) dx-/c ^ c-^^cfvkdxv'' (A = 0,1, 2, 3). 



Considering first the projection of p.2.7p in the direction of m, we remark that one may use (j3.1.6p 
and (|3.1.10[) to conclude that for solutions, 7/^9^6^'' = is equivalent to equation (|4.1.ip . 

We now project (|3.2.7p onto the orthogonal complement of u, which, with the aid of (|3.1.8p . gives 
the three equations n{,9^9'"' = 0, j = 1, 2, 3 : 

(3.2.21) = n^a^e^"- = Ili[R, + c-^P] {e^u'')d,,{e^u'') + (ni9»^(n,^ - ^2^) 

- ni [Rc + c-^P] {e^u^d^ie^u") + {Uid''^){Rc - 3c-^P). 

After making the substitutions (|3.2.12p . ()3.2.13p . (|3.2.14p . and (|3.2.15p . and using relation (|3.2.20p . 
it follows that for solutions, (|3.2.2ip is equivalent to (|4.1.3p . 

We also introduce the nameless quantity Qc and make use of (|3.1.10p . (|3.1.12p . (|3.1.14p . p.l.l5p . 
(|3.1.16p . ([33^ . (|3X5)) . and (|3.2.12p to express it in the following form: 

(3.2.22) Q/=^n^ = = 0.(.,p, $), 
where 

(3.2.23) £lc{v,p, $) &l{r,,p)e^^/^' [^c{^,p) + c'^ = &l{v,P)[Rc + c-'P]- 
Then we use the chain rule together with p.l.6p . (|4.1.ip . and (|3.2.22p to derive 

(3.2.24) e^u^'a^F + Qcd^^ie^u'^) = (4F - 3Qc)e^u''a^</), 

which we may use in place of (|3.1.6p . Upon making the substitutions (|3.2.4p . p.2.5p . (|3.2.12p . 
(I3.2.13P . and p.2.14p . and using the relation (|3.2.20p . it follows that for solutions, (|3.2.24p is 
equivalent to (|4.1.2p . 

4. The Formal Limit c ^ oo of the EN^ System 

For convenience, in this section we list the final form of the EN^ system as derived in sections 
13.11 and 13.21 We also take the formal limit c ^ c» to arrive at the EP^ system and introduce the 
equations of variation (EOV|;). 

4.1. A Recap of the EN^ System. 

The EN^ system is given by 

(4.1.1) dtri + v''dkV ^0 

(4.1.2) dtP + v'^dkP + Qcdkv'' + c-^jcfQcVk{dtv'' + v^'Oav'') 

= (4P - 3Qc) [c-^dt'S> + c-^v''dk<^>] 

(4.1.3) (7c)'(i?c + c-^P)[dtv^ + v'^dkv^ + c-\-fc)^v^Vk{dtv'' + v^'dav'')] + d,P 

+ C-\lcfv^{dtP + v'^dkP) = {ic-^P - i?c)(aj* + (7,)-2^;^■[c-2a,$ + C-^v'^dk'^]) 

(4.1.4) - c-'^dl^ + A$ - = AttG{Rc - 3c-^P), 
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where j = 1, 2, 3, 
(4.1.5) 

7c = 7c (v) 



(c2 - |v|2)l/2 



(4.1.6) 



i?/=^'e4*/c^5H,(r;,p) 



(4.1.7) 



^='Oe(ry,p,$) ^=^S^(r?,p)e4*/^^[5H,(,7,p) + c-2p] = ( ^{v-.p) e^^/'^' [m^p) + c-'p] 



(4.1.8) 

P^^e^^/^'p, 



c denotes the speed of Hght, 6^1], p), which is defined in (|3.1.16p . is the speed of sound, and the 
functions and 6c derive from a c— indexed equation of state as discussed in Section l3.ll The 
variables 77,p,v = {v^ ,v'^ ,v^), and $ denote the entropy density, pressure, (Newtonian) velocity, 
and cosmological- Nordstrom potential respectively. Section [6] contains a detailed discussion of the 
c-dependence of the ENJ: System. 

4.2. The EPk System as a Formal Limit. 

Taking the formal limit c ^ cx) in the ENJ: system gives the Euler-Poisson system with a cosmological 
constant: 



(4.2.1) 
(4.2.2) 

g22I) 

(4.2.3) 
(4.2.4) 
where 
(4.2.5) 

(4.2.6) 



dtv + v^dkTl = 
dtp + v^dkP + Qoodkv'' = 



(j = 1,2,3) 



^$ = 47rGi?„ 



Qc 



3^oo(??,p) 



Mp) = 61 



•j {V,P)^oo{V;P) 



dp 



■(»7,P) 



9^oo(??,p), 



9loo('7,p) and &l^{r],p) are the limits as c ^ oo of DXciv^P) and 62(?7,p) respectively (see (|6.3.ip . 
(I6.3.2p . and (|6.3.3p '). and the quantity Rao is the Newtonian mass density. Since equations (|3.1.16p 
and (I6.3.3P imply that dDlociVi p)/9p = 6^{r],p), it then follows with the aid of the chain rule that 
for solutions, equations (|4.2.2p and (|4.2.2[ ) are equivalent. We refer to the solution variable $ 
from equation (I4.2.4P as the cosmological-Newtonian potential. 

This system of equations is discussed in |10| . in which, under an isothermal equation of state {p — 
c^p^ where the constant Cg denotes the speed of sound), Kiessling derives the Jeans dispersion relation 
that arises from linearizing (|4.2.2f ). (|4.2.3p . (|4.2.4p about a static state in which the background mass 
density p is non-zero, followed by taking the limit k — > 0. 

It is a standard result that the solution to (|4.2.4p is given by 



(4.2.7) 



G 



-k|s — s'l 



s - s' 



[fnoo(r/(t,s'),p(t,s')) -5Hoo(?7,p)] dV, 



where the constants ^oo^V, and p, which are the values of rj, and p respectively in a constant 
background state, are discussed in Section [S] The boundary conditions leading to this solution are 
that $(t, •) — l>oo vanishes at oo, and we view $(t, s) as a (not necessarily small) perturbation of the 
constant potential $oo- 

Remark 4.1. Consider the kernel /C(s) = — Ge~''l^l/|s| appearing in (|4.2.7p . An easy computation 
gives that /C(s), V'^^^/C(s) G L-'^(M''). Therefore, a basic result from harmonic analysis (Young's 
inequality) implies that the map f ^ K, * f, where * denotes convolution, is a bounded linear 
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maiJ^ from L'^{M.^) to H^iSJ^). From this fact and Remark IB. 21 (alternatively consult Lemma [6-1 p . 
it follows that $(t, ■) e H^+^(R^) whenever fj{t, ■),p{t, •) e Hj^ (R^), {M.^) respectively. By then 
applying Lemma IA-41 we can further conclude that •) G H^+'^{R^) whenever T]{t, ■),p{t, ■) £ 

5. The Equations of Variation (EOV^) 

The EOV^ are formed by linearizing the EN)^ system (EP^ system in the case c = oo) around a 
background solution (BGS) V of the form V = {rj, P,v^ , ■ ■ ■ , $2, ^3)- Given such a V and inhomo- 
geneous terms f,g, h^^\ h^'^\ h'^^\ I , we define the EOVJ^ by 



(5.0.8) dtri + v'^dkV - / 

(5.0.9) dtp + v'^dkP + Qcdkv'' + c-^^,fQ,Vk{dtv'' + v^dav'') = g 

(5.0.10) (7c)'(i?c + C-^P) [dtv' + v'^dkV^ + C-2(^j2~j~ .fe ^ yag^^k^^ 

+djP + c-\%fv^{dtP + v^dkP) = h^^'^ 
(5.0.11) -c~^d^^ + A<i>-K^^ = l, 



where 7c c/(c^ — |wp)^/^, i?c "= e^*/^^5Hc(?7,p); ^tc. The unknowns are the components of 
W '= {i),P,v'^,v^,v^) and <i>. 

Remark 5.1. We place parentheses around the superscripts of the inhomogeneous terms h'-^^ in 
order to emphasize that we are merely labeling them, and that in general, we do not associate any 
transformation properties to them under changes of coordinates. 

Let us now provide a few remarks on our notation. We find it useful to analyze both the dependent 
variable p and the dependent variable P when discussing solutions to (|4.1.ip - (|4.1.4p . Therefore, 
we will make use of all four of the following arrays: 



(5.0.12) W={'n,P,v\v^,v^) 

(5.0.13) V = (r,, P, v\v^,v\ $, 9^$, 82^, 83^) 

(5.0.14) W={tj,p,v\v^,v^) 

(5.0.15) V = v\v^, $, 82^, 93$), 



where P e''*/'^ p. When discussing a BGS V (77, P, u^, ■ • ■ , <i>2, $3) that defines the coefficients 
of the unknowns in the EOV^, we also use notation similar to that used in (|5.0.12p - (|5.0.15p . 
including V {rj,p,v^, ■ ■ ■ ,83^), W =^ {rj,P,v^,v^,v^), where p =^ g-4*/c ^^^^ When c — 00, 
we may also refer to W {rj,p,v^ ,v^) as the BGS, since in this case, the left-hand sides of 
(|5.0.8p - (|5.0.1ip do not depend on and furthermore, W = W. Additionally, we may refer to the 
unknowns in the EOV^ as IV {fi,p,v^ ,v'^ ,v^) when c = 00; in this article, <i> will always vanish 
at infinity, and in the case c = 00, rather than considering $ to be an "unknown," we assume that 
the solution variable $ has been constructed via the convolution $ = A^] * /, where the kernel A^(s) 
is defined in Remark 1 4. 1[ and / is the right-hand side of (|5.0.1ip . 

We frequently adopt standard PDE matrix notation. For example, we may write (|4.1.ip - (|4.1.3p 

as 

(5.0.16) c^^(W,$)a^W = b, 

where each c-^'^ { ) is a 5 x 5 matrix with entries that are functions of W and <&, while b = 
{f,g,--- , /i^^-*) is the 5-component column array on the right-hand side of (|4.1.ip - (|4.1.3p . 

Remark 5.2. We emphasize that throughout this article, we operate under the convention that 
the c-A'^i') are functions of the BGS variables W, $. We therefore write "cj4''('W, $)," as opposed to 
writing ^A^iW,^)." 



'Our proof breaks down at this point in the case k = 0. 
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It is instructive to see the form of the c-A'^i'), — 0, 1,2,3, for we will soon concern ourselves 



with their large— c asymptotic behavior. Abbreviating ac 



he? (Re 



2p^ 



(5.0.17) c^°(W,$) 



/I 











\ 





1 

















o(l-3) 









a,(l + /3p)) 


fl(2,3) 


Vo 


Qc0^ 


fl(3,l) 




ac(l + /3PV 



(5.0.18) 



=^"(W) 



/I 











^ 





1 

















Roo 

















Roo 





^0 











^oo / 



(5.0.19) cA\W,<^) 



(5.0.20) 







































l + /3p) 


acV 















pi'''' 






) 






Vo 










a 


,1;1(1+/3PV 








/v^ 








\ 










Qoo 












) = 






1 Rocv^ 





















RooV^ 















i^O 


i? 




V 





etc. 



6. On the c-Dependence of the EN^ System 

In addition to appearing directly as the term c~^, the constant c appears in equations (|4.1.ip - 
(I4T4)) through four terms: 1) P = e'**/'=V, 2) 7^ = c/{c'^ - |v|^)^/^ 3) Rc = e'**/'='jnc(?7,p), and 
4) Qc = Sc(^;-P)s"'*^^' [^c(f7,p) + c~^p]. Because we want to recover the EP^ system in the large 
c limit, the first obvious requirement we have is that the function fHc('7,p) has a limit 9loo('7,p) as 
c — > 00. For mathematical reasons, we will demand convergence in the norm | • Iat+i^it at a rate of order 
c^^, where £ is a compact subset of M+ x R+ that depends on the initial data; see (|6.3.ip and (|6.3.2p . 
Although a construction of £ from the initial data is described in detail in Section 18.21 let us now 
provide a preliminary description that is sufficient for our current purposes: for given initial data, we 
will prove the existence of compact sets O2, 02, [—a, a]^, if =^ (92 x [— a^]^, =^ D2 x [—a, a]^, and 
a time interval [0, T] so that for all large c, the (c— dependent) solutionO V (V) to the EN^ system 
launched by the initial data exist on [0, T] xR^ and satisfy W([0, T] xM^) c 62, W([0, T] xR3)^c D2, 
y([0,T] X R3) c K, and V([0,T] x R^) ^ ^ gee Section lO for a detailed description of 62 and 
D2, and (|10.2.25p . (|10.2.26p for the construction of K and K 

The set £ from above, then, is the projection of D2 onto the first two axes (recall definition 
(15.0. 14p ). Intuitively, we would like the aforementioned four functions of the state-space variables 
to converge to p, l,i?oo, and Qoo respectively when their domains are restricted to an appropriate 
compact subset. In this section, we will develop and then assume hypotheses on the c— indexed 
equation of state that will allow us to prove useful versions of these kinds of convergence results. 



^Recall the notation 15.0.1211 - 15.0.1511 which defines the arrays W, V, "W, and V respectively. 
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6.1. Functions with c Independent Properties: The Definitions. 

The main technical difficulty that we must confront is ensuring that the Sobolev estimates pro- 
vided by the propositions appearing in Appendix [B] can be made independently of all large c. By 
examining these propositions, one could anticipate that this amounts to analyzing the norms (see 
definition (|2.3.4p ) of various c— indexed families of functions arising in the family of EN^ systems. 
We therefore introduce here some machinery that will allow us to easily discuss uniform-in-c esti- 
mates. Following this, we use this machinery to prove some preliminary lemmas that will be used 
in the proofs of Theorem [2] and Theorem lU which are the two main theorems of this article. Before 
proceeding, we refer the reader to the notation defined in (j2.4.ip . which will be used frequently in 
the discussion that follows. 

Definition 6.1. We define TiP {c^\ D;qi, - ■ ■ , qn) to be the ring consisting of all c— indexed families of 
functions of qi , • • • , such that for all large c, £ (S), and such that the following estimate 
holds: 

(6.1.1) <c^--C(S)). 

We emphasize that the constant C(S)) is allowed to depend on the family and the set £>, but 
within a given family and on a fixed set, C(S)) must be independent of all large c. 

Remark 6.1. At the beginning of Section Wl^ we explain why it is sometimes useful to shift the 
point of view as to what are the arguments of '^c- Therefore, is perhaps best thought of as a 
"c— indexed expression" rather than a "c— indexed family of functions." 

Definition 6.2. Assume j > 2. Then {c'';'D; qi, - ■ ■ , qn) denotes the sub-ring contained in 
7?.^ (c'^'; £); gi, • • ■ , g„) consisting of all such c— indexed functions ^dqi, ■ ■ ■ , Qn) such that 

(6.1.2) Uciqir-- ,qn)\\H. <c'' -CiD-.WqiW^, , • • • , ) 

holds for all constant arrays (gi, ^2, • • • , 9n) G S and all array- valued functions (gi,--- ,g„) G 
H^^R^) X ••• X i?|„(M^) having {(gi(s), g2(s), • • • ,g„(s)) | s e M^} ^ rpj^g constant C{T)) is 
allowed to depend on the family and the set S), but it can depend on the gi, • • • , g,i only through 
their H^. normj0. We remark that the spaces are defined by (|2.3.1[) . 

Remark 6.2. This definition is highly motivated by (jB.33[) of Appendix [Bl 

Remark 6.3. We also emphasize that in our applications below, the qi and qi may themselves depend 
on the parameter c, even though we do not always explicitly indicate this dependence. Typically, the 
qi will be quantities related to solutions of the EN^ system, and the qi will be equal to the components 
of either ([8T2| . (jS.l.lOp . or (|8. 1.111) . perhaps scaled by a power of c. 

Notation. If G V {c^] S; gi, • • • , g„), then we sometimes write 

(6.1.3) ^c{qi.--- ,g„) ^©^(c^Digi,--- ,g„). 

Remark 6.4. We employ the following abuse of notation throughout this article: if we have written 
"5'c('?i,''' i^n) £ (c*"; 3D; 91, • • • ,qn)" then we are indicating that in addition to the c'^— type 
bound on 5'c given in (|6.1.ip , that the functions qi have the properties stated in Definition 16.21 
i.e., we use the notation in quotations to also communicate that S) is compact and convex, that 
there are constants qi, which will be clear from context, such that qi e H^.(M.^), that the image set 
{(gi(s), g2(s), • • • ,g„(s)) | s e R^} is contained in D, and that {qi,q2,--- ,g„) G T). We employ a 
similar abuse of notation in writing "5^c('Zi, • • ■ , qn) G Z^(c'^; S; gi, • ■ • , qn)" 

Remark 6.5. In the notation TZ{- • ■ ),!(• ■ • ), and 0-'(- ■ • ), we often omit the argument D. In this 
case, it is understood that there is an implied set 2) that is to be inferred from context; frequently 
D is to be inferred from L°° estimates on the qi that follow from Sobolev embedding. Also, we omit 
the argument c*^ when fc = 0. Furthermore, we have chosen to omit dependence on the constants 
iji since, as will be explained at the beginning of Section 16.31 their definitions will be clear from 
context. We will occasionally omit additional arguments when the context is clear. 

"'^^Technically speaking, the || • arc not norms in general, since HOH^fj ^. = oo unless qi = 0. This is not a 

problem because in this article, we only study the || ■ "norm" of functions q{s) that by design feature < oo. 
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6.2. Functions with c Independent Properties: Useful Lemmas. The foUowing three lem- 
mas provide the core structure for analyzing the Sobolev norms of terms appearing in the EN^ 
system. They are especially useful for keeping track of powers of c. Their proofs are based on the 
Sobolev-Moser estimates that are stated as propositions in Appendix IbI 

Lemma 6-1. If j > 2 and ^dqi, ■ ■ ■ , Qn) G TZ^{c''-D; gi, • • • , q,i), then 
^c{qi,--- ,qn)~dc{qi,--- ,q,i) ^V{c'';'Z);qi, - ■ ■ ,g„). 



Proof. We emphasize that the conclusion of Lemma 16-11 is exactly the statement that 
||g'c(gi,--- , 9n)-5'c(gi, • • ■ ,9n)||ffj l^cf" -CiWqiWfji Jknll/fj ). Its proof foUows from definitions 
O and 121 and from (|R33|. ' " □ 

Lemma 6-2. Suppose that j > 2, g'c G W {c^^^D^qi, ■ ■ ■ ,g„), C5c e (c^^^ ; 2); gi, • • • ,g„), 
andiDc eXJ(c'^^;D;gi,-- - ,g„). Then {^c ■ &c){qu- ■ ■ , g„) G 7^J (c'=i+'^^ ; S; gi, • • • ,g„) 
and {dc-^c){qi,--- , g„) G (c'=i+'=3; D; g^, . . . 

Proof. Lemma [6-21 follows from the product rule for derivatives and (|B.30p . □ 

Remark 6.6. Lemma [6-21 shows that for fc < 0, TZ^{c'';'D; qi, - ■ ■ , qn) is a ring i.e., it is closed under 
products. We frequently use this property in this article without explicitly mentioning it. 

Remark 6.7. Lemma [6-21 can easily be used to show that ii G ^"'(c°; 2); 9i, • ' ' j and for all 
large c, doesn't vanish on £), then 1/dc G "/^■'(c"; T);qi, - ■ ■ , qn). 

Remark 6.8. Lemma 16-21 shows that if dc{qi,''' ,qn) = 0, then G X-' (c*^ ; S ; 51 , • • • In 
particular, if g = 0, then any polynomial (of strictly positive degree) in q is an element of V (q). 

Remark 6.9. Lemma [6-21 shows in particular that for k < 0, 2^c'^;T);qi, ■ ■ ■ ,(?„) is an ideal in 

7^J(^;gl,■•• ,qn). 

Remark 6.10. li k < and there is a fixed function 5^00 G 'TV{D\ qi, - ■ ■ , qn) such that — doo G 
W {c'^;^;qi, ■ ■ ■ ,qn), then it follows that |t?c|j,s ^ ItJooIj.s + 1, so that the family is uniformly 
bounded in the norm | • for all large c. A similar remark using the || • norm applies if 
doo G I-'(S); • ■ • ,qn) and ^c — ^oo & P {c'';T);qi, ■ ■ ■ ,qn)- We often make use of these observations 
in this article without explicitly mentioning it. 

Lemma 6-3. Suppose that j > 3, ki + k2 ~ ko, and that Sc{qi, • • ■ , qji) G W {c'^°; Di; qi, - ■ ■ , qn)- 
Assume further that for 1 < i < n, G C'^([0, T], Hq.) D C^([0, T], Hq. ) and that for all large c, 
0*^2 {dtqi, ■■■ , dtqn) ([0, T] x M?) C Ss. Then on [0, T], we have that 



dt 



dciqi,- ■ ■ ,qn) i(c'=i;S)i X D2;gi, • • • ,g„,c'=2(9tgi, ■ • ■ ,c'==^9t9„ 



Proof. Lemma 16-31 follows from the chain rule. Lemma 16-21 and Remark 16.81 We emphasize that 
the constant term associated to c^'^dtqi is 0, so that on the right-hand side of the definition (|6.1.2p 
of ■ • ), we are measuring c'^'^dtqi in the H^^^{M.^) norm. □ 

Corollary 6-4. Let da be a first- order spatial differential operator. Suppose that j > 3, ki + k2 = ko, 
and that Sc{qi, ' ' ' , qn) G TV {c^°] Si; (71, • • • , qn)- Assume that for all large c, 
c^^{daqi,--- ,9ag„)([0,r] xM?) c£)2. Then on [0,T], we have that 

da(dc{qi, ■ ■ ■ ,9n)) G p-'^{c''^;Di X D2]qi, ■ ■ ■ ,qn,c''^daqi, ■ ■ ■ ^cJ^'-daqn)- 

Proof. The proof of Corollary 16-41 is virtually identical to the proof of Lemma [6-31 □ 

6.3. Application to the ENJj. System. 

We will now apply these lemmas to the ENJ: system. Let us first make a few remarks about our 
use of the norms || • that appear on the right-hand side of (|6.1.2p and the constant term qi 

associated to qi. For the remainder of this article, it is to be understood that the constant term 
associated to c'^V is c*''Vc, that the constant term associated to c^V is c^Vci and the constant 
term associated to both DV and DV is 0, where Vc and Vc are defined in l|8.1.10p and (|8.1.1ip 
respectively. Furthermore, the constant term associated to c^'^V is understood to be c~'^Vc, and 
so forth. In other words, when estimating c'^V using a th order Sobolev norm, it is understood 
that we are using the norm || • , and similarly for the other state-space arrays. The relationship 

c'=V 

between the arrays V and V is always understood to be the one implied by (|5.0.13p and (|5.0.15p . 
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We furthermore emphasize that V (or V) wiU represent a solution array to the EN^ system, and 
therefore will implicitly depend on c through the c— dependent initial data Vc and through the c 
dependence of the ENJ^ system itself. The fact that the constant arrays Vc and Vc depend on the 
parameter c does not pose any difficulty. For as we shall see, Vc is contained in the fixed compact 
set A for all large c, and is contained in the fixed compact set K for all large c, where the sets 
K and K were introduced at the beginning of Section [SI Therefore, when we require L°° estimates 
of the constants Vc and Vc, the bounds can be made independently of all large c. 

In addition to the above remarks, we add that we will have available a-priori estimates that 
guarantee that V e C°([0, T], iJ^JnCHP, T],H^-^) nC^{[0, T], H^'^) for a fixed integei0iV > 4 
on our time interval [0, T] of interest, which are hypotheses that are relevant for Lemma 16-31 and 
Corollarv l6-4l Our a-priori estimates will also ensure that all of the relevant quantities are contained 
in an appropriate fixed compact convex set, so that the "hypotheses on the q.;" described in Definition 
I6.2l will always be satisfied. Consequently, we will often omit the dependence of the running constants 
C(- • ■ ) (see Section [2?5|) on such sets. The relevant a-priori estimates ("Induction Hypotheses") are 
described in detail in Section riO.3.11 

Let us now provide a clarifying example and also note that as we change settings, it is sometimes 
useful to shift the point of view as to what are the arguments of a family ^c{- ■ ■ )■ For example, 
consider the expression ^c == c~^9t$, where $ is a solution variable in the EN^ system depending 
on c through the initial data Vc and through the c— dependence of the system itself. If it is known 
that c""'^ ||(3f$||^3 is uniformly bounded by L for all large c, then we have that ^c G 2^{c~'^; c~^9t$) 
since ||c~^9t$||//3 < c~^L. If it also turns out that ||9t$||//3 is uniformly bounded for all large 
c, then have that ^c G I^ic~^;dt^). If both estimates are true, then we indicate this by writing 
5c el^{c-^;c^^dt^)nl^{c-^;dt<^) or 5c = 0^{c-^;c-^dt<^)r]0^{c-^;dt'i>). These kinds of estimates 
will enter into our continuous induction argument in Scction ri0.21 in which we will first prove a bound 
for c~^5f<i>, and then use it to obtain a bound for dt^; see (|10.2.21l) and (|10.2.23p . 

Remark 6.11. For simplicity, we are not always optimal in our estimates. 

The following four lemmas, which provide an analysis of the c— dependence of the terms appearing 
in the EN^ system, will be used heavily in Section 110.31 which contains most of our technical 
estimates. Before providing the lemmas, we first restate our hypotheses on the equation of state 
using our new notation. 

Hypotheses on the c Dependence of the Equation of State. 

(6.3.1) mdv.p), ^^ooiv,?) e n''+\€;v,p) 

(6.3.2) d\civ,P)-'^oc{v,P) e 7^^+^(c-^;e:;r,,p), 

where the set £ was introduced at the beginning of Section [5] and is described in detail in Section 
18.21 We emphasize that our construction of £ will depend only on the initial data, and not on the 
equation of state (although it is clearly not the case that (|6.3.ip and (|6.3.2p are satisfied by all 
c— indexed equations of state). As a simple consequence of p.l.l6p . (|6.3.ip . and (|6.3.2p . we have 
that 

(6.3.3) eliv^p) - eUv.P) e n'^ic-^; C; ri,p). 

We also assume that 9^oo(?7,p) and &'^{r],p) are "physical" as defined in Section [23 i.e., we assume 
in particular that whenever rj^p > 0, we have < 9loo(?7,p) and < 6^{r/,p). 

Hypothesis (|6.3.ip ensures that the terms appearing in the EN^ and EP^ systems are sufficiently 
differcntiablc functions of V, thus enabling us to apply the Sobolev-Moser type inequalities appearing 
in Appendix [Bl It is strong enough to imply Theorem [1] and Theorem [31 Hypothesis (|6.3.2p is used 
in our proof of Theorem [2] and Theorem [D Although a weakened version of Hypothesis (|6.3.2p is 
sufficient to prove a convergence theorem, we do not pursue this matter here since we are not striving 
for optimal results. 

Lemma 6-5. With the functions 7c, Rc, Roa, Qc, Qoo, W, and W of the state-space variables defined 
in (jiXS]) . (jiX^ . (gXij), ((iX7|) . (jiXS)) . (|5.0.12p . and (|5.0.14p respectively, we have form ^ 0, 1,2 



'The relevance of Af > 4 is described in Section |8] 
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and = t, 1, 2, 3 the following estimates for all large c, including c = oo : 

(6.3.4) (^j2_^^^iv+i(^-2.^) 

(6.3.5) e^*/^' - 1 e 7^^+l(c'"-2. c-™$) 

(6.3.6) Rc-Roo = e4*/^'5H,(77,p) - 5Hoo(??,p) e 7^^+l(c'"-2; r,,p, c-"<i>) 

(6.3.7) Qc - Qoo = Qc{v,P, - noo(?7,P) e 7^'^(c"-'; r?,p, c~"$) 

(6.3.8) W - W e 7^^(c'"-2;P,c""$) 

(6.3.9) W e 7^^(W,c-"$) 

(6.3.10) a^w - a^w e I^-\c"'-^;P, d^P, 0-"*$, c-'"^^*) 

(6.3.11) e Z^-i(W, V(i)W, c-'"$, c-'"a^$). 

Proof. (|6.3.4p . and (|6.3.5|) are easy Taylor estimates. (|6.3.6p follows from Lemma|631 (|6.3.ip . (j6.3.2p . 
and (jO^l) . (jgXT)) then follows from p.l.l6p . p.2.23p . g^H), Lemma [El and (jgX^ . 

Since F - p (1 - e-'**/=')P, ([QIS]) follows from ([eXS]) . Lemma [631 and that the fact that W 
and W differ only in that the second component of W is p, while the second component of W is P. 
(I6.3.9P is a simple consequence of (16.3.81) . (|6.3.10p follows from (I6.3.8p . Lemma [6-31 and Corollary 
(|6.3.1ip then follows easily from (|6.3.10p . □ 



Lemma 6-6. IfO<j<N,k<0, and e W{c'';W), then for m = 0, 1, 2, we have that 

(6.3.12) 5c e 7^^(c'^■;W,c-"$). 

Lemma [6-6l follows easily from expressing W in terms of W and c^™$ via (|6.3.9p and applying the 
chain rule. Let us re-phrase the content of Lemma 16-61 in order to be clear: assuming the quantity 
5c, when expressed in terms of the state-space variables TV, has the regularity/c'"'— boundedness 
properties defined by the class of functions TIP [c^] TV), then upon re-expressing the same quantity 5c 

in terms of the variables W, $, ^i.e. ^cCW) = 5c('W(W, $))^ , it has the regularity/c'"'— boundedness 

properties defined by the class of functions 7^^ (c'^'; W, c^™$). 

Lemma 6-7. Let c.A^CW, $), v — 0,1,2,3, denote the matrix-valued functions ofW and $ intro- 
duced in Section\^ Let the c— dependent relationship between TV and W, $ be defined by (j5.0.12p 
and ()5.0.14p . Then for all large c including c = oo, and for m = 0, 1, 2, we have that 

(6.3.13) oo.A''(TV), (oo.A"(TV))"^ e 7^^(TV) n7^^(W,c-™$) 

(6.3.14) c.A''(TV,$), (cyi°(TV,$))"^ e 7^^(TV,c~™$)n7^^(W,c-™$) 

(6.3.15) c^''(TV, $) - oo.A'^(TV) e 7^^(c™-2; TV, c"™*) n 7^^(c™-2. w, c"™*) 

(6.3.16) {cA°iW,^)y^ - (o„yi°(TV))"^ en'^ic""'^;W,c-"'^)r\n^{c'''~'^;W,c-"'<P). 

Proof (|6.3.13P - (|6.3.16P follow from (|5.G.17p - (|5.G.20p . Remark[6Jl Lemma[E2l Lemma[E5l Lemma 
16-61 the determinant-adjoint formula for the inverse of a matrix, and the hypotheses (|6.3.ip . (I6.3.2P 
on the equation of state. □ 

Lemma 6-8. Let ^^{WM^^<i>) (o, 0, -$Koo(ry,p)ai$, -$Koo(?7,p)a2$, -5Hoo(r/,p)a3$)) denote 
the right-hand side of (|4.2.ip . (|4.2.2p . (|4.2.3p . and let Q3c(TV, $, £>$) rfenoie t/ie right-hand side 
()4.1.ip - (|4.1.3p . Let the c^ dependent relationship between TV and W, $ &e defined by (|5.0.12p and 
(|5.0.14p . T/ieTT. for all large c including c = oo, and for m — 0,1, 2, we have that 

(6.3.17) «Boo(TV,v(i)$) eX^(TV,V(i)$) nZ^(W,c-™$, V(i)$), 

(6.3.18) »c(TV, D$) e I^(TV, 0"™$, V^^^^, c-"5t$) n I^(V^, c""$, v'^^^, c-™at$), 
and furthermore, 

(6.3.19) !Be(TV, $, i^$) = «B„o(TV, V^^^*) + 5c, 
w/iere 

(6.3.20) 5c e X^(c"-2; W, c"'"*, £-"£»$). 
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Proof. (16.3. 17P - (|6.3.20|) all follow from combining the fact that *Bc(Wc, 0, 0, 0, 0) = with 
Remark 16.81 Lemma 16-21 Lemma 16-51 and Lemma 16-61 □ 



Example 6.1. As an enlightening example, we discuss the non-relativistic limit of the polytropic 
equation of state, that is, an equation of state of the form p — mgc^n + ^^^^n^, where 7 > 1. 

Let us assume that Ac.Aoo e $n^+Hni(G:);?7), that A^o > on ni(e:), and that A^ - A^o e 
7?.^+^(c^^; ni((r); ry), where ni(£) is the projection of the set £ introduced at the beginning of 
Section [6] onto the first axis. Some omitted calculations show that Hypotheses 16.3.1) and I6.3!2l then 

') P 



hold, and that 






(6.3.21) 


Rc 


Al'\r,) 




(6.3.22) 


Qc 




(6.3.23) 






(6.3.24) 







c2(7-l) 



In the isentropic case ry(t, s) = 77, (|6.3.23p can be rewritten in the familiar form p = C • (Roo)'^ , 
where C is a constant. 

7. Energy Currents 

In this section we provide energy currents and discuss their two key properties: 1) for a fixed c, 
they are positive definite in the variations W when contracted against certain covectors, and 2) their 
divergence is lower order in the variations. In Section [8.31 '^6 will see that the positivity property is 
uniform for all large c. A general framework for the construction of energy currents for hyperbolic 
systems derivable from a Lagrangian is developed in [6J. The role of energy currents is to replace 
the energy principle available for symmetric hyperbolic systems by providing integral identities, or 
more generally, integral inequalities, that enable one to control Sobolev norms of solutions^ to the 
EOV^. This technique will be used in our proofs of Lemma [10-141 and Theorem |4l 

7.1. The Definition of an Energy Current. 

Given a variation W : — > and a V : — + 'M}^ as defined in Section |5l we define 

the energy current to be the vectorfield ^^^^3 with components ^'^'>3 , ^'^^S' , j = 1,2,3, in the global 
rectangular coordinate system given by 

(7.1.1) 

f,^ + C + 2c-^^l{v,v'^)P + t{Rc + c-'P) [v.v'^ + c-'j'Av.v'f] 

(^)3' ^:!=!Wrf -f ^P^ + 2 [v^ + c-^^lvW] P + jURc + cr'Py [vkv" + c-'j^A^kv'^f] ■ 
In the case c = oo, we define for j = 1, 2, 3 : 
(7.1.2) ^3 =z 71 ~\- — \- K^VkV 

Qoo 

dcf ^ 2vip + R^vHkV^. 

Qoo 

We note that formally, linic^oo ^'^^d = a fact that will be rigorously justified in Section [5751 

The energy current (|7.1.ip is very closely related to the energy current J introduced in [18j , where 
the following changes have been made. First, we have dropped the terms from J corresponding to 
the variations of the potential $ and its derivatives, for we will bound these terms in a Sobolev norm 
using a separate argument. Second, the expression for ^^'3 is constructed using the velocity state- 
space variable v (|3.2.1ip and variations v, as opposed to the variables e'^u^ and variations 



"'^^As we shall see, the energy currents ^'^^3 do not control the variations <i? or Di; these terms are controlled 
through a separate argument that uses the lemmas in Appendix 1X1 
^^Recall that we also refer to W as the BGS when c = oo. 
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il^ that appear in the expression for j. Finally, we emphasize that the formula for '"^^3 applies in 
a rectangular coordinate system with = t, whereas in the formula for J'^ provided in [18j . the 
rectangular coordinate system is such that — ct, even though c was set equal to unity in |18| . 

Remark 7.1. Viewed as a quadratic form, is manifestly positive definite in the variations W 

if p > 0, for by our fundamental assumptions on the equation of state, p > i?oo > and 

Qoo > 0. 

7.2. The Positive Definiteness of for ^ e 

As discussed in detail in [TB], for ^ belonging to a certain subset of T*M, the quadratic fornj^ 
f^/'^''3^(W, W) is positive definite in W if P > 0. The energy current j from [T8| has the property 
that ^^J^ is a positive definite quadratic form in V for ^ belonging to the interior of the positive 
component of the light cone in T*M. In contrast, since the energy current '■'^^3 from (|7.1.ip does 
not contain terms involving the variations of the potential <i>, '^^^^^-'^''(W, W) is positive definite in 
{•f],P,v) for ^ belonging to 1^*'^ , the interior of the positive component of the sound cone at x, 
which is larger than the light cone. Expressed in coordinates, this statement reads 



(7.2.1) ^^('=)/(W,W) > if W > 0, P>0, and^eZ|*+, 
where 

(7.2.2) Tr+ = {C G T*M I (h-^y^C^C. < and Co > 0}, 

and is the reciprocal acoustical metri^^ with components that read (for j, fc = 1, 2, 3) 

(7.2.3) h"'^-c-^~i^,)'[e-^i?,,p}-c-^] 

(7.2.4) = = -(7,)2 [e-^i?i,p) - c-^]v^ 

(7.2.5) h^'' = 5'^ - (7c)' [e^^{v,p) - c-^]v^v^ 



in the global rectangular coordinate system. Recall that the function ©c is defined in | 3. 1.141 

We now further discuss the reciprocal acoustical metric. The characteristic subselm of the trun- 
cate43 EOV^ (IS.O.Sp - (jS.O.lOp . which is a subset of T*M., the cotangent space at x, is the union 
of several sheets. The inner sheet is the sound cone at x, which is expressed in coordinates as 
{C G T*A4 I {h~'^)^'^ C,^C,i, = 0}. It follows from the general construction of energy currents as pre- 
sented in that ^^^'^•'^'^(W, W) is positive definite whenever ^ belongs to the interior of the positive 
component of the inner sheet, which is exactly the condition expressed in (|7.2.1|) . This fact allows 
us to use the form ■f^^'-^'S (W, W) to estimate the norms of the variations W, provided that we 
estimate the BGS V. 

As an alternative justification of the fact (|7.2.ip . we remark that ^^'3 has the same form as an 
energy current which is shown to have the sound cone positivity property by Christodoulou [7], 
except that instead of using Christodoulou's rectangular coordinate system on Ai featuring = t, 
we are using here a rectangular coordinate system with x'^ — ct. In addition, we have expressed ^^^3 
using the Newtonian velocity v and its associated variations v, as opposed to Christodoulou's use 
of the four velocity and its variations. 

Remark 7.2. Because limc^oo &c'^iv,'P) — &^{'n,'P) > 0, it follows that for all large c, the covector 
with coordinates (1, 0, 0, 0) is an element of I^*^ ■ Therefore, ('^)3*'(W, W) is positive definite for all 
large c. 



"'^'^We write"^^('^^3''(W, W)" to emphasize the point of view that ^(j'^^'S^ is a quadratic form in W. 
^'^The reciprocal acoustical metric was introduced using dimensionless variables in [7j. 

^^|18| contains a detailed discussion of the notion of the characteristic subset in the context of the ENJ^"^ system. 
■^■^By "truncated EOV^" we mean the system that results upon deleting the variable # and the equation 115.0. lll l 
that it satisfies. 
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7.3. The Divergence of the Energy Current. 

As described in [TSj, if the variations W are solutions of the EOV); (I5.0.8P - (|5.0.10p . then we can 

compute dfj,(^^'^^d'^^ and use the equations (jS.O.Sp - (|5.0.10[) for substitution to eliminate the termfl 
containing the derivatives of W : 



(7.3.1) 



p2 



+ 2c-2a)4 



Rc 



c-'P 



2f]J + 2^ + 2vjh^^\ 

Qc 



We observe here that in the case c — oo, (|7.3.ip reduces to the more palatable expression 



(7.3.2) 



dt 



Qc 



f + {dt{Rco) + dj{RaoV^)}vkv'' 



2rif + 2^ + 2vjh^^l 

Qoo 



8. Assumptions on the Initial Data and the Uniform-in-c Positivity of Energy 

Currents 

In this section we describe a class of initial data for which our energy methods allow us to 
rigorously take the limit c ^ c» in the EN^ system. The Cauchy surface we consider is 
{{t,s) €M\t = 0}. 

8.1. An ff^ Perturbation of a Quiet Fluid. 

Initial data for the EP^ system are denoted by 



.1.1) 



Voo(s) = (?7,p,w^^^^>^$oo,*o,^'l,*2,*3), 



where ^1/0(8) = dt^(0,s) and "^fj — 9j$oo(s). We assume that Voo is an H perturbation of the 
constant state Voo, where 

(8.1.2) Voo = {fj,p, 0, 0, 0, $00, 0, 0, 0, 0), 

fj,p are positive constants, and the constant $00 is the unique solution to 

(8.1.3) K^^^+'inGm^ifj,p)=0. 

The constraint ()8.1.3p must be satisfied in order for equation ()4.2.4p to be satisfied by Voo- By an 
perturbation, we mean that 

(8.1.4) I|Woo||h« <oo, 

where we use the notation Woo and Wqo to refer to the first 5 components of Voo and Voo respec- 
tively. We emphasize that a further positivity restriction on the initial data p and 77 is introduced in 
Section \8.2l and that throughout this article, N is a fixed integer satisfying 



.1.5) 



TV > 4. 



Remark 8.1. We require A^ > 4 so that Corollarv lB-31 and Remark IB. II can be applied to conclude 
that Of I e C°{[0,T],H'^~^{R^)), where I is defined in (|10.2.10p : this is a necessary hypothesis for 
Proposition I A-5l which we use in our proof of Theorem [2l 



^■^Showing this via a calculation is an arduous task. The lower-order divergence property is a generic feature of 
energy currents constructed in the manner described in 6 , but we require its explicit form in order to analyze its c— 
dependence. 
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Although we refer to $00 and ^J^^, ly — 0,1, 2, 3, as "data," in the EP^ system, these 5 quantities 
are determined by fi,p,v^ through the equations (|4.2.2[ '). (|4.2.4p . and (|8.1.3p . together with 
vanishing conditions at infinity on ~ ^00 and ^0 ■ 

(8.1.6) A$oo - '*'($oo - $00) = 47rG[inoo(?7,p) - ^oo(?7,p)] 

(8.1.7) A*o - «'*o = -47rG5t|t=o(^oo(r/,p)) = -\T^Gdk[^^{f],v)v'') , 

where the integral kernel from (|4.2.7p can be used to compute $00 — ^00 and ^q. We will nevertheless 
refer to the array Voo as the "data" for the EP^ system. 

Remark 8.2. Remark SH] imphes that $00 e and *oo e 

We now construct data for the EN^ system from Voo. Depending on which set of state-space 
variables we are working with, we denote the data for the ENJ^ system by 

(8.1.8) V, = (77,p, v\v', i„ ^-n, ^2, ^-3) 

(8.1.9) or V, (77, e4*=/^'p, v\v' ,v\ 4„ *o, ^1, ^-2, ^s), 

where unlike in the EP^ case, <lc, ^ii ^2, and are data in the sense that the system is under- 
determined if they are not prescribed. We have chosen the data T],'p, ,iP' , xr" , ^q, ^'2, ^3 for the 
EN^ system to be the same as the data for the EP„ system, but for technical reasons described below 
and indicated in (|8.1.12p and (|8.1.14p . our requirement that there exists a constant background state 
typically constrains the datum <i>c so that it differs from <i>oo by a small constant that vanishes as 

C — > CXD. 

As in the EP^ system, we assume that Vc is an perturbation of the constant state of the 
form (depending on which collection of state-space variables we are working with) 

(8.1.10) V,= (^,p, 0,0,0, $,,0,0, 0,0) 

(8.1.11) or Ve = [j], Pc, 0, 0, 0, i-„ 0, 0, 0, 0) 

where fj and p are the same constants appearing in Voo , '^c is the unique solution to 

(8.1.12) K^l., + 47rGe''""'*<= [9^c(?y,p) - ic^^p] = 0, 

and Pc "= e^'^ ^^"p. The constraint (|8.1.12p must be satisfied in order for equation (|4.1.4p to be 
satisfied by p, fj, and Although the background potential $c for the EN^ system is not in general 
equal to the background potential $00 for the EP^ system, it follows from the hypotheses (|6.3.ip 
and (|6.3.2p on the c-dependence of that 

(8.1.13) lim$e = $oo. 

c — >oo 

We now define the initial datum $2 appearing in the arrays (|8.1.8p and (|8.1.9p by 

(8.1.14) ic = ioo - ^00 + $c, 

which ensures that the deviation of ic from the background potential $c matches the deviation of 
$00 from the background potential $00. We denote the first 5 components of V^,, Vc, Vc, and Vc 
by Wc, WcWc, and Wc respectively. 

Remark 8.3. We could weaken the hypotheses by allowing the initial data for the EN^ system to 
deviate from the initial data for the EP^ system by an perturbation that decays to rapidly 
enough as c — > 00. For simplicity, we will not pursue this analysis here. 

8.2. The Sets 0,02,02,€,K, and J?. 

In order to avoid studying the free boundary problem, and in order to avoid singularities in the 
energy currents (|7.1.ip and (|7.1.2p . we assume that the initial pressure, energy density, and speed 
of sound are uniformly bounded from below by a positive constant. According to our assumptions 
p.l.lip on the equation of state, to achieve this uniform bound, it is sufficient to make the following 
further assumption on the initial data: that 'VVoo(K'^) is contained in a compact subset of the 
following open subset O of the state-space , the admissible subset of truncated state-space, defined 

by 
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(8.2.1) 0={W^{r],p,v\v^,v^) e M'"^ | 77 > 0,p > O}. 

Therefore, we assume that WocCM"^) C Oi and Woo G Oi, where Oi is a precompact open set 
with Oi (s O, and " (s " means that "the closure is compact and contained in the interior of." We 
then fix convex precompact open subsets O2 and D2 with Oi (s ©2 ^2 ^ O, and define £ to 
be the projection of D2 onto the first two axes, where O2 denotes the closure of D2. We assume 
that with this definitior0 of £, hypotheses (|6.3.ip and (|6.3.2p are satisfied by the equation of state. 
Consequently, property (|8.1.13p shows that for all large c including c = 00, Wc(K'^) (s O2 and 
Wc e Int(02); also note that for aU c, W^ = Woo = Woo. 

We now address the variables ($, 9t$, (92*&, c?3*&) ■ In Section [TOl we will use energy esti- 
mates to prove the existence of an interval [0,T] and a cube of the form [— a,a]^ such that for 
all large c including c = cx), we have ($, 92^", 93$) ([0, T] x R^) C [— a,a]^. Furthermore, 

it will follow from the discussion in Section [10] that for all large c including c = 00, we have 
($c, *o, ^1, ^2, ^3)(IR^) <£ Int([-a, a]^). The compact convex set K, then, as given in (|10.2.25p 
below, will be defined to be O2 x [— a, a]^. It follows from the above discussion that for all large 
c including c = cx), we have Vc(M^) d Int(_ft') and Vc G lnt{K). Our goal is to show that the 
solution Vc to (|4.1.ip - (|4.1.8p launched by the initial data Vc exists on a time interval [0, T] that 
is independent of (all large) c and remains in K. 

We now discuss the simple construction of ^ : based on the above construction, it follows from 
definitions (|5.0.12p - (|5.0.15p that for all large c including c — 00, we have Y G K =^ W G D2. 
As given in (|10.2.26p . we will then define the compact convcjj3 set Rhy M.'^ O2 x [—a, a]^, so that 
for all large c including c = 00, we also have V G A' => V G As in the previous discussion, it 
follows that for all large c including c = 00, we have Vc(K3) (e Int(j?) and Vc G Int(il). 

8.3. The Uniform-in-c Positive Definiteness of '■'^^3 ■ As mentioned at the beginning of Section 
[71 we will use the quantity (OIIl^ to control ||W(i)||^jv, where ^^'3 is an energy current for the 

variation W with coefficients defined by a BGS V. Since we seek estimates that are uniform in c, it 
is important that '^^■'3 is uniformly positive definite in W independent of both the BGS V and all 
large c. Let us now formulate this precisely as a lemma. 

Lemma 8-1. Let be the energy current ()7.1.ip /or the variation W defined by the BGS V. 
Assume that W(i,s) G 02, where O2 is defined in Section \8.2[ and furthermore assume that 
1^ Z. Then there exists a constant Cq^ ^ with < Cq^ z < 1 such that 

(8.3.1) Co,,z\^? < (^)a"(w,w) < c^ljwi' 

holds for all large c including c = 00. 

Proof. It is sufficient prove inequality (|8.3.ip when |W| = 1 since it is invariant under any re-scaling 
of W. Let W,V be the arrays related to the arrays W,V as defined in ()5.0.12p - ()5.0.15|) . Our 
assumptions imply the existence of a compact set £> depending only on O2 and Z such that for all 
large c, V{t,s) G D. 

Recall that is defined in (|7.1.2p and that is manifestly positive definite in the varia- 

tion0 W if p > 0. If we view as a function of (W, W), then by uniform continuity, there is 

a constant < C{D) < 1 such that C(2))|W|2 < < C(S))~^ |W|2 holds on the compact set 

{|W| = 1} X 115(2)), where 115(2)) is the projection of T) onto the first five axes. Furthermore, if 
we also view ^^^3 as a function of (W, V), then by Lemma [6-21 Lemma [6-51 (17.1.11) . and (|7.1.2p we 

^"^Note that our construction of O2 depends on 'Woo(R'^), but is independent of the equation of state and inde- 
pendent of c. 

•^^Proposition IB^ requires the convexity of K and .8, and the estimate IIB.331 1 also requires that Vc K, V £ A 
In practice, K and A can be chosen to be cubes. 

■^'^To be consistent the notation used in formula 117. 1.21 1 , it would be "more correct" to use the symbol W to denote 
the variations appearing as arguments in '°°'3(-,-). However, for the purposes of this proof, there is no harm in 
identifying W = W since in this context, these placeholder variables merely represent the arguments of when 
viewed as a quadratic form in the variations. 
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have that (^'j = '■°°^d + S'c ■ where € 7^^(c-2; D; V). (|8XT|l now easily follows: C^^.z 

can be any positive number that is strictly smaller than C(£'). □ 

9. Smoothing the Initial Data 

For technical reasons, we need to smooth the initial data. Without smoothing, the terms on the 
right-hand sides of (|10.2.4p - (|10.2.6p involving the derivatives of the initial data could be unbounded 
in the norm. To begin, we fix a Fricdrichs moUifier i.e., x G C^(M'^), 
supp(x) C {s| |s| < l},x > 0, and / x d^s = 1. For e > 0, we set Xe(s) =' e^'^x(f)- We smooth the 
first 5 components Woo of the data Voo defined in (jS.l.ip with Xe, defining Xe'^oo G C°° by 

(9.0.2) x.Woo(s)'='/ Xe{s-s')W^.{s')d's'. 

Note that we do not smooth the data (l>c, ^o) G H^^'^ xH^^^ because by Remark lST^ and definition 
(|8.1.14p . they already have sufficient regularity. 

The following property of such a mollification is well known: 

(9.0.3) lim ||x.Woo-Woo||ff« =0. 

e— »0+ 

We will choose below an eq > 0. Once chosen, we define 



(9.0.4) (°) W ((°)^, (°)v) Xe„ Woo 

(9.0.5) (°) ((0)^, e**^/^' . v) , 

where <i>c is defined in (|8.1.14p . By Sobolev embedding, the assumptions on the initial data Wc, 
which are the first 5 components of the data Vc defined in (|8.1.9p . by Lemma [6-21 by (|6.3.5p . and 
by the mollification property (|9.0.3p . 3{Ai > A eo > 0} such that 

(9.0.6) for all large c, ||W - '^^^'W^Wh'^ < Ai =^ W G O2 

(9.0.7) ||(o)W,-W,||^« <C0^,^.^, 

where O2 is defined in Section [8T2l and Cq^ z is the constant from (|8.3.ip . Here, Z is a fixed constant 
that will serve as an upper bound for ||$(t)||icxj on a certain time interval, where $ will be a solution 
variable to the EN^ system. We explain this fixed value of given in expression (|10.3.5p below, in 
detail in Section flO. 31 Note that according to this reasoning, Ai = Ki{02', Z). 

Remark 9.1. Because these quantities enter into our Sobolev estimates below, it is an important 
fact that ||(°^Vc||„iv+i and ||(°)Vc|1loc are uniformly bounded for aU large c. By (I8.1.13p . (|8.1.14p . 

definition (j9.0.5p . and Sobolev embedding, to obtain these uniform bounds, we only need to show 
that II 6^*'=/'^ • '-^^p ||^f"+i is uniformly bounded for all large c. This fact follows from Lemma 

16-11 Lemma [6-21 and (|6.3.5p . Such a uniform bound is used, for example, in the estimate (|10.3.45p . 

10. Uniform-in-Time Local Existence for EN^ 

In this section we prove our first important theorem, namely that there is a uniform time interval 
[0,T] on which solutions to the ENj: system having the initial data Vc exist, as long as c is large 
enough. We emphasize that throughout this article, we assume that > 4 and that the equation 
of state satisfies the hypotheses (16.3.11) . (|6.3.2p . 

10.1. Local Existence and Uniqueness for EN^ Revisited. 

Let us first recall the following local existence result proved in [TS] , in which it was not yet shown 
that the time interval of existence can be chosen independently of all large c. 

Theorem 1. (EN^ Local Existence Revisited) Let Vc(s) he initial data (|8.1.8p for the EN'^^ 
system (j4.1.ip - (|4.1.8p that are subject to the conditions described in Section O Assume that 
the equation of state is "physical" as described in Section \ 3.1\ . Then for all large (finite) c, there 
exists a Tc > such that (|4.1.ip - (|4.1.8p has a unique classical solution V S C^([0,Tc] x M.^) of 
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the form V — {■r],P,v^,v'^,v^,^, 81^,81^,82^,83^) with V(0,s) = Vc(s). The solution satisfies 
V([0, Tc] X K'^) C K, where the (c— independent) compact set K is defined in equation (110. 2. 25^ . 
Furthermore, V G ^([0, TJ, H C^P, T,], i/^-^) f\C'^{[0,T^],H^-'^) and 
$ e CO([0,T,],i/|'+i)nCi([0,rd,i7|'JnC2([0,Te]^i/|'-i)nC3([0,Te],V|'-2), «;/iere i/ie constants 
Vc and are defined by (jS.l.lip and (j8.1.12p respectively. 

Remark 10.1. Although they are not explicitly proved in the facts that V e C^([0, TJ x M^) 
and that V is twice differentiable in t as a map from [0, Tc] to -ff-^"^ follow from our assumption 
that N > A (i.e. the degree of differentiability of the solution is TV — 2). Also, by Corollary IB-3[ we 
have that p e C°([0, T^], F^) n Ci([0, TJ, iJ^"^) n C2([0, TJ, ffp^-^), since p = Pe-^^/c'. 

The proof of the claim that Tc can be chosen such that V([0, Tc] x M^) c K is based on the fact 
Vc(K'^) (E Int(iir), which follows from the remaining discussion in this section, together with the 
continuity result from the theorem and Sobolev embedding. 

Remark 10.2. The case c = 00 is discussed separately in Theorem [31 

Remark 10.3. The local existence theorem in [18^ was proved using the relativistic state-space 

variables L/" =^ e'^u'^. However, the form of the Newtonian change of variables made in sections l3.ll 
and 13.21 together with Corollary IB-3[ allow us to conclude Sobolev regularity in one set of variables 
if the same regularity is known in the other set of variables. 

The following corollary, which slightly extends the lifespan of the solution and also allows us to 
conclude stronger regularity properties from weaker regularity properties, will soon be used in our 
proof of Proposition 110-21 

Corollary 10-1. Assume that 5/2 < N' < N and that V(t,s) is a solution to the EN^. system 
having the regularity property V G L°"{[{),T\,H^ ) fl C"'^([0,T] x M^). Then there exists an e > 

such that ( with T+ = T + e) 

(10.1.1) V e C^{[0,T+] X R^)nc'^i[o,T+],H^jnc\[o,T+],H^-^)nc^{[o,T+],H^-^). 

Proof. We apply Thcorcm[l]to conclude^ that for each T' G [0, T], there exists an e > 0, depending 

on T', and a solution V to the EN^ system such that V G C'^{[T' - e,T' + e] x R^) 

nCO([r ~€,r + e],H^) n C\[T' ~€,r + e],H^^^) n C2([r' ~e,T' + e],H^-^) and such that 

V(r') = V(r'). Furthermore, the uniqueness argument from [18 , which is based on local energy 
estimates, can be easily modified to show that solutions to the ENJ: system are unique in the class 
C^{[T' — e,T' + e] x M'^). Therefore V = V on their common slab of spacetime existence. Corollary 
[TO thus follows. □ 

In addition to Theorem [TJ our proof of Theorem [2| also requires an additional key ingredient, 
namely a continuation principle for Sobolev norm-bounded solutions: 

Proposition 10-2. (Continuation Principle) Let Vc(s) be initial data ()8.1.8p for the EN'^ 
system ()4.1.1|1 - ()4.1.8|1 that are subject to the conditions described in Section [3 and let T > 0. 
Assume thatV (E C^{[Q,T)xR^)r\C°{[0,T),H^ )nC^{[0,T), H^^'^) is the unique classical solution 

existing on [Q,T) launched by Vc(s). Let O be the admissible subset of truncated state-space defined 
in ()8.2.ip . and let lis : R^^ — > R^ be the projection onto the first 5 axes. Assume that there are 
constants Mi,M2 > 0, a compact set K C M^" with I{q{K) d O, and a set U d Lnt(K) such that 
the following three estimates hold for any T' G [0, T) : 

(1) III V |||ff«^,T'< Ml 

(2) III 8tV \\\hn-i j^,< M2 

(3) V([0,T'] X M^j c [/. 

Then there exists an e > such that ( with T^ = T + e) 

V G ( [0, r+] X n c° ( [0, r+] , i/^j n ( [0, T+] , h^- ^)nc^{[o,T+], h^-^) 

(10.1.2) and V([0,T+] x R^) C K. 



•^^Theorem [T] can be easily modified to obtain a solution that exists both "forward" and "backward" in time. 
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Remark 10.4. Hypothesis (2) is redundant; it can be deduced from hypothesis (1) by using the 
equations to solve for dtV and then applying (|B.30p . 

Proof. We will first show that there exists a V* e i?^ (M^) such that 

(10.1.3) lim ||V(T„)-V*||h«-i =0 

n — ^oo 

holds for any sequence {Tn} of time values converging to T from below. 

If {Tn} is such a sequence, then hypothesis (2) implies that ||V(r,) - V(Tfe)||H"-i < M2\Tj-Tk\. 
By the completeness of H^~^, there exists a V* G H^~^ such that (I10.1.3|) holds, and it is easy 
to check that V* does not depend on the sequence {Tn}. By hypothesis (1), we also have that 
{V(T„)} converges weakly in to V* and that ||V*||^jy < Mi. We now fix a number N' with 

5/2 < N' < N. By Proposition [B^ we have that lim„^oo1|V(r„) - V*||^„' = 0. Consequently, 
if we define V(r) V*, it follows that V £ L°°{[0,T], H^J n C"{[0,T], H^'j. Using the fact 
that N' > 5/2, together with the embedding of (K^) into appropriate Holder spaces, it can be 
shown that V £ C°{[0,T], H^') =^ V, V^^^V G C°{[0,T] x R^); i.e., we can continuously extend 
V,V(i)V to the slab [0,r] x R^. To conclude that V e C^{[0,T] x R^), we wih show that dtV 
extends continuously to [0,T] x R^. To this end, we use the ENJ: equations to solve for 9tV : 

(10.1.4) 9tV = 5(V,V(i)V), 

where d G C"^. Since V, V^^^ V G C°([0, T] x R^), the right-hand side of p0.1.4p has been shown to 
extend continuously so that it is an element of C°([0,T] x R'^). Furthermore, since 
V G C^([0,T) X K^), it follows from elementary analysis that dt'V exists classically on [0,T] x R^ 
and that dtV G C°([0, T] x R^), thus implying that V G C^dO, T] x R^). The additional conclusions 
in (|10.1.2p now follow from Corollary 110-11 and continuity. □ 



Remark 10.5. Proposition 1 1 0-21 shows that if the solution V blows up at time T, then either 
ImiT'^T III V \\\[{N rpi— oo, \\iB.T''\T III dtV |||ff"-i.T'— oo, or V(T',R'^) escapes every compact 
subset of O X i?5 as T' ] T, where O is defined in ((8XT|) . 

Remark 10.6. Although the main theorems in this article require that N > A, Corollarv 110-11 and 
Proposition [TO^ are also valid for = 3, except that the conclusion V G C'^{[0,T + e] x R'^) must 
be replaced with V G C^{[Q,T + e] x M^), and the conclusion V G C'^{[0,T + e],H^-^) does not 
hold. 

10.2. The Uniform-in-Time Local Existence Theorem. 

We now state and prove the uniform time of existence theorem. 

Theorem 2. (Uniform Time of Existence) Let Voo denote initial data (IS.l.ip for the EP^ 
system (j4.2.ip - (j4.2.6p that are subject to the conditions described in Section O Let Vc denote 
the corresponding initial data (j8.1.9p for the EN"^ .system (|4.1.ip - (|4.1.8p constructed from Voo «s 
described in Sectionl^ and let (°^Wc denote the smoothing (j9.0.5p of the first 5 components ofVc o,s 
described in Section\^ Assume that the c— indexed equation of state satisfies the hypotheses (j6.3.ip 
and (|6.3.2p and is "physical" as described in sections \3.1\ and 1 6. SI and let K be the fixed compact 
subset o/ M^*^ defined in (jl0.2.25p . Then there exist cq > and T > 0, with T not depending on 
c, such that for c > cq, Vc launches a unique classical solution V to (j4.1.ip - (|4.1.8p that exists 
on the slab [0,T] x R^ and that has the properties V(0,s) = Vc(s) and V([0,r] x M^) c K. The 
solution is of the form V = (r/, P, , v^, v^, $, 9t$, dicf), 92$, ^3$) and has the regularity properties 
VeC\[0,T] xR^)nC°{[0,T],H^^{R^))nC\[0,T],H^-\R'^)) and 

$ G C3([0,T]xR3)nC"([0,T],i?|'+i(M3))nCi([0,T],i/£(M3))nC72([0,r],i/£-i(R3)), where the 
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constants Vc and $c o.'i^e defined by (|8.1.1ip and (|8.1.12p respectively. Furthermore, with 
p pg-40/c ^ ffiQj-g exist constants Ai, A2, Li, L2, is, L4 > such that 



(10.2.1a) III W-(*"WJ||^iv,i.< Ai 

(10.2.1b) III $ - $c |||//"+i_T< A2 

(10.2.1c) III i9tW |||^iv-i_-r< 

(10.2. Id) III llljjiv L2 

(10.2. le) III 9(^?7 lll^jv-2 III i9^p lll/i-iv-2 r< L3 

(10.2. If) c"^ III 9(2$ |||//«-i^T< ^4- 



10.2.1. Outline of the structure of the proof of Theorem\^ We prove Theorem [5] via the method of 
continuous induction ("bootstrapping"). After defining the constants Ai,A2,L2, and £4, we make 
the assumptions (|lG.3.ip - (|lG.3.4p . These assumptions hold at t = and therefore, by Theorem [l] 
there exists an interval r G [0, T^) on which the solution exists and on which the assumptions hold. 
Based on these estimates, we use a collection of technical lemmas derived from energy estimates to 
conclude that the bounds (|10.2.17p - (|10.2.23p hold for r € [0, Tc). It is important that the constants 
appearing on the right-hand sides of (|10.2.17p - (|10.2.23p do not depend on c, if c is large enough. 
We can therefore apply Proposition ll0-2l to conclude that for all large c, the solution can be extended 
to a uniform interval [0, T]. The closing of the induction argument is largely due to the fact that the 
source term for the Klein-Gordon equation satisfied by which is the right-hand side of (|4.1.4p . 
"depends on <f> only through c^'^^." 

10.2.2. Proof of Theorem\^ We begin our detailed proof of Theorem [2] by making a few remarks 
about the running constants C( - ■ • ) that will appear in our argument. For the sake of appearances, 
we frequently suppress the dependence of C(- ■ • ) on iV, k, and the sets ©2, D2, and ^. We indicate 
dependence of the running constants on the initial data ||'^")Wc||j:^iv+i, |j<i>c||j:^iv+i , and ||^o||//« by 

writing C{id). By Remark l9. 11 any constant C{id) can be chosen to be independent of all large c. 

We now introduce some notation that will be used throughout the proof, and also in the following 
section, where we have placed the proofs of the technical lemmas. Let V denote the local in time 
solution to the ENJ^ system (|4.1.ip - (|4.1.8p launched by the initial data Vc as furnished by Theorem 
[TJ With W denoting the first 5 components of V, we suggestively define 

(10.2.2) W(<, s) W(t, s) - (°)Wc(s) 

(10.2.3) $ = $-$c, 

where 4^ is defined in (|8.1.14p and (°)Wc(s) is defined in (pUS]) with the help of (|10.3.5p . We 
remark that this choice of (°^Wc(s) is explained in more detail below. 

It follows from the fact that W is a solution to (|4.1.ip - (|4.1.3p that W is a solution to the 
EOV^ (inm - (I5.O.IOP defined by the BGS V with initial data W(0,s) = Wc(s) - (°)Wc(s). 
The inhomogeneous terms in the EOV^ satisfied by W are given by b = (/, g, • • • , ft.^"^-*), where for 
J -1,2,3 

(10.2.4) / = -i;'-afep77] 

(10.2.5) 5 = (4P - 2,Qc)[dt{c-^<^) + v'' dk{cr^ <^)] - v'^ dk[e^^^'''" ■ 

(10.2.6) /i^^^ = (3c-2p-i?c)(9j3' + (7c)"^w^'[at(c-2$) +„'=aj^(c-2$)]) 

In order to show that the hypotheses of Proposition 110-21 are satisfied, we will need to estimate 
in L^. Therefore, we study the equation that 9sW satisfies: for < |d| < we differentiate 
the EOVJj defined by the BGS V with inhomogeneous terms b to which W is a solution, obtaining 
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that (9qW satisfieo 

(10.2.7) c^^(W, $)9;.(95W) = ba, 
where (suppressing the dependence of the c-^^i-) on W and <&) 

(10.2.8) hs = e^"55 ((c^°)-'b) + ks 
and 

(10.2.9) ks = [(cyL")-^yi'afc(9sW) - 95((c^°)~'c^'5fcW)] . 

Thus, each S^W is a solution the EOV^ defined by the same BGS V with inhomogeneous terms 
hs- Furthermore, $ is a solution to the EOVJ: equation (|5.0.1ip with ^{t = 0) = 0, and the 
inhomogeneous term I on the right-hand side of (|5.0.1ip is 

(10.2.10) / = {k^ - A)i, + 4ttG{Rc - ic-^P). 

We will return to these facts in Section 110.31 where we will use them in the proofs of some technical 
lemmas. 

As an intermediate step in our proof of (|10.2.1aP - (jl0.2.1f|) . we will prove the following weaker 
version of (jl0.2.1dp : 

(|10.2.1d[ ) \\\ dt^ \\\hn L'^- 

We now define the constants Ai,A2, Lj, and L4. We will then use a variety of energy estimates 
to define Li, L2, and L3 in terms of these four constants and to show that p0.2.1ap - (|10.2.1f|) are 
satisfied if T is small enough. First, to motivate our definitions of L'2t L-^, and A2, see inequalities 
(IA.4P and (|A.6p of Lemma IX^ and inequality (|A.20p of CoroUarv lA-3[ and let Co(k) denote the 
constant that appears throughout the lemma and its corollary. By a non-optimal application of 
Lemma ri 0-61 we have that 

(10.2.11) Co{K){c-'\\ifo\\H- + WmU--^) < 1/2 

(10.2.12) Co(k)(c||/(0)||^«-i + ||(A - k^)4/o - dtm\\„N-2) < 1. 

Note also the trivial (and not optimal) estimate {Co{k))^c^^\\'^q\\]jn < 1/4 == (A2)^/4. With these 
considerations in mind, we define 

(10.2.13) A2''=^'l 

(10.2.14) L'2 = l 

(10.2.15) Li = l. 

To define Ai, we first define Z = Z{id; A2) to be the constant appearing in (|10.3.5p . Using this 
value of Z, which we emphasize depends only on A2 and the initial data Woo for the EP^ system, 
we then define Ai so that (|9.0.6p and (|9.0.7p hold. Note that it is exactly at this step in the proof 
that the smoothing ^"^Wc, which is defined in (|9.0.5p . of the initial data Wc, which are the first 5 
components of (|8.1.9p . is fixed. 

We find it illuminating display the dependence of other constants that will appear below on 
Ai, A2, L2, L4. Therefore, we continue to refer to (|10.2.13p - (|10.2.15p by the symbols A2, ^2; L4 
respectively, even though they are equal to 1. 

We now carry out the continuous induction in detail. Let J"™"^ be the maximal time for which 
the solution exists and satisfies the estimates (110.2. lap . (110.2. Ibl) . (|10.2.1d[ ). and (|10.2.1fp : i.e., 

J.n^ax '^^^ ^up {t\Y g C"^ {[0 , T], H^^iR^)) H C\[0, T], H^^\R^)) H C\[0, T], H^-^M."^)) , 

(10.2.16) and ()10.2.1ap . (|10.2.1bp . (|10.2.1df ). and (|10.2.1fp holdj. 

Note that the set we are taking the sup of necessarily contains positive values of T since for all large 
c, the relevant bounds are satisfied at T = 0, and therefore by Theorem [H also for short times. 



^^Recall the convention stated in Remark l5.2l 
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Lemmas [T^Tl [TO^ [T^Tl \T0AT\ [T^9l and inequalities (|10.3.30p and (|10.3.29p of Lemma [T^TOl 
supply the following estimates which are valid for < r < T"''"'^ : 

(10.2.17) III W |||H«,r< [Ai/2 + T . C(Ai, A2, L,,L'^)] ■ exp(T • C(Ai, A2, Li, L'^)) 

(10.2.18) III dtW |||^«,,<Li(Ai,A2,L'2) 

(10.2.19) III d^rj |||h«-2,.,III d^p\\\H---^r<L3i^iA2,Li,L'^,U) 

(10.2.20) III $ |||^„+i^,< ^^+C7(Ai,A2,L2)+T2.C(Ai,A2,Li,L'2,L3,i4) 

(10.2.21) |||5t$ |||H«,r<4/2 + r-C(Ai,A2,Li,4) 

(10.2.22) III 5(2$ |||ffiv-i,^<i4/2 + r-C(Ai,A2,Li, 4,^3,^4) 

(10.2.23) 1119*$ \\\H".r< i2(Ai, A2, ii, 4)72 + T • C(Ai, A2, Li, 4, L3, L4). 

We apply the following sequence of reasoning to interpret the above inequalities: first Li in 
(110.2. ISp is determined through the known constants Ai,A2, and L'2. Then L3 in p0.2.19p is de- 
termined through the known constants Ai,A2,Li,L2, and L4. Then L2 in (|10.2.23p is determined 
through Ai, A2, Li, and £'2- Finally, the remaining constants C(- • • ) in p0.2.17p - (I10.2.22p are all 
determined through Ai, A2, Li, Lj, L3, L4. 

By Sobolev embedding and (|8.1.13p . there exists a cube [— a, a]^ (depending on the initial data, 
Ai, and L2) such that for all large c, || $ ||//iv< Ai, || 9*$ ||hn< -^2 =^ 

(10.2.24) (^$, ^2$, <93$, 9*$) ([0, T] x R^) C [-a, a]^ 

independent of whether or not $ or 9f$ are solutions to an equation. Motivated by these consider- 
ations, we define both for use now and use later in the article the following compact sets: 

(10.2.25) K = 62x[-a,a]'^ 

(10.2.26) K = D2x[-a,a]^. 

Here, O2 and D2 are the sets defined in Section [8?2l 

We now choose T so that when < r < T, it algebraically follows that the right-hand sides of 
(|10.2.17p and (|10.2.20p - (|10.2.23p are strictly less than Ai, (A2)^ ^4, and L2 respectively. Note 
that T may be chosen independently of (all large) c. We now show that T""*^ < T is impossible. 

Assume that T™"^ < T. Then observe that the right-hand sides of (|10.2.17p and (|10.2.20p - 
(|10.2.23p are strictly less than Ai, (A2)^, L2, L4, and L2 respectively when r — T^°-^. Therefore, by 
the construction of the set K described above, by (|9.0.6p . and by Sobolev embedding, we conclude 
that for all large c, Vc([0, T™"^) x M."^) is contained in the interior of K. Consequently, we may apply 
Proposition 110-21 to extend the solution in time beyond T™"^, thus contradicting the definition of 
ymax^ Note that this argument also shows that V([0, T] x R^) c K. This completes the proof of The- 
orem [5] □ 

10.3. The Technical Lemmas. 

We now state and prove the technical lemmas quoted in the proof of Theorem [2] We will re- 
quire some auxiliary lemmas along the way. Throughout this section, we assume the hypotheses 
and notation used in our proof of Theorem [21 i.e., V denotes the solution, W denotes its first 5 
components, the relationship between W and W is given by (|5.0.12p and (|5.0.14p . W and $ are 
defined in ()10.2.2p and (|10.2.3p respectively, I is defined in ()10.2.10p . and so forth. In this section, 
we typically suppress the dependence of the running constants C(- • • ) on N, k,02, O2, K, and ^. 

10.3.1. The induction hypotheses. We assume that r e [0,r™''^). By the definition (|10.2.16p of 
r™°^, we have the following bounds (|10.3.2p - (|10.3.4p (we will comment on the bound (|10.3.ip in 
a moment), where Ai,L2, and L4 are defined in (|10.2.13p - (|10.2.15p respectively: 

(10.3.1) III W-(")W, |||^«,,< Ai 

(10.3.2) III -S-lc |||ff«+i,,< A2 

(10.3.3) c'' \\\dt^\\\H«.r<L'2 

(10.3.4) III dt<i> \ \\h»---,t< U- 
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We note the foUowing easy consequence of (|8.1.14[) and (|10.3.2p : 

(HnXl) III $ - $c \\\H^,r <lll $ - *c |||ff«^r + 1 1 1 " ^^c ||U«.r< A2 + C{id} = C{ld; A2). 

It then follows from (|8.1.13p . (|8.1.14p . (110.3. 2[ ). and Sobolev embedding that 

(10.3.5) III $ |||lco_^< Z(id; A2). 

Using the value of Z in (110. 3. Sp . which depends only on the data Woo for the EP^ system and the 
known constant A2, we are able to choose a constant Ai > such that ()9.0.6p and (j9.0.7p hold. As 
discussed in sections [51 and [ro.2.2l such a choice of Ai also involves fixing the smoothing ("^W of 
Woo, which then defines ^^^Wc via equation (|9.0.5p . We emphasize that it is this choice of ^^^Wc 
and Ai that appear imphcitly in (|10.2.17p . in (|10.2.1ap . and in (|10.3.ip . 
Then by ((5117)) and (|10.3.ip . we have that 

(HnXB) 

IIIW-W, |||^«,, < |||W-(°)W, |||ff«,, + III ("'W.-W, |||^«,, + III We-We|||^^«,, 

< Ai+C{id;Ai) + C{id) C{id; Ai). 
Furthermore, by Lemma [6-11 (|6.3.9p . and (|10.3.1f ). we have that 

(10.3.6) III W- Wc \\\H«,r< C{id;Ai,A2). 

Note also that (|9.0.6p . (|10.3.ip . and the definition of D2 given in Section [57^ together imply that 
for all large c, we have that W([0,T™'^^) x R^) C 62 and W([0, T^"""') x R^) ^ 

In our discussion below, we will refer to (|10.3.ip - (|10.3.6p p0.3.1[ ). and (|10.3.2[ ) as the induction 
hypotheses. Sobolev embedding and the induction hypotheses, which for all large c are satisfied at 
T = 0, together imply that W, V'^^W, W, V^^^ W, $, V^^^*, c'^dt^, c~^df^ are each contained in a 
compact, convex set (depending only on the initial data, Ai, A2,L2j ^^.d L4) on [0,T™''^) x M^. As 
stated in Remark 16.51 we will make use of this fact without explicitly mentioning it every time. 

10.3.2. Proofs of the technical lemmas. 

Lemma 10-3. Consider the quantity I defined in (I10.2.10p . Then for m = 0, 1, 2, we have 



(10.3.7) 


{4nG)-H = y\oo{v,P) - ^Hoo(?7,p) +;?c, 


(10.3.8) 


{47TG)-^dtl ^ dt{mo,{v,P)) + 


(10.3.9) 


{47TG)-'d^l = df{D\oo{v,P))+^c, 


where 




(10.3.10) 


;?cex'^(c"-2;^,p,c-™$) 


(10.3.11) 


G5c e X^^i(c"-2; c"™*, dtV, dtP, c~"'dt<^) 


(10.3.12) 





Proof. It follows from the discussion in Section [5] that 
(10.3.13) 

Therefore, (|10.3.7p + (|10.3.10p follows from Lemma O Lemma |S31 and Lemma O (|10.3.8I) + 
(110.3. lip and p0.3.9p + p0.3.12p then follow from Lemma [O □ 

Lemma 10-4. 

(10.3.14) III dtW \\\HN-^,r, III dtW \\\H«-Kr< G {td; A, , A2 , L'^) =^Li(jd; Ai, A2, 4). 

Proof. By using the ENJ: equations (|4.1.ip - (|4.1.3p to solve for 9tW and applying Lemma 16-21 
(|6.3.10p in the cases v = 1,2,3, Lemma Lemma and Lemma we have that 

(10.3.15) 

dtW = (e^°(W, $))"' [ - ,A''{W, <f)dkW + ^.(W, $, D$)] 

= (oo^°(W))"'[-oo^'=(W)5fcW + 25oo(W,v(i)$)] +0^-i(c-2;W,V«W,<i>,V«$) 
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The bound for ||| (9tW |||ijiv-i ^ now follows from Lemma [6-2[ (16.3. 13p . (|6.3.17[) . the induction 
hypotheses, (|10.3.15[) . and the definition of 0^^^(- • • ). The bound for ||| dt'W |||^iv-i ^ then follows 
from the bound for 1 1 1 9tW 1 1 \h^-'^,t: l|6.3.10p in the case ly — t,m = 1, and the induction hypotheses. 
We remark that we have written the "intersection term" on the right-hand side of (|10.3.15p in a 
form that will be useful in our proofs of Lemma [10-61 and Lemma [10-71 the "c~^" decay is used in 
Lemma 110-61 while the "dependence on c~^Z?$" is used in Lemma 110-71 Similar comments apply 
to Corollary 110-51 and equation (|10.3.18p below. □ 

The following indispensable corollary shows that for large c, the EN^ system can be written as a 
small perturbation of the EP^ system. 

Corollary 10-5. (EN^ « EP«, for Large c) 
(10.3.16) 

dtW (oo^°(w))"' [ - ooA\w)dkW + ^oo{w, v(i)$)] + o^-i(c-2; w, v(i)w, v(i)$) 

Proof. Recall that dfW and dt'W differ only in that the second component of dfW is dtP, while the 
second component of c^jW is dtP- Therefore, it follows trivially from (|10.3.15p that (|10.3.16p holds 
for all the components of (JtW except for the second component dtP- 

To handle the component dtp, we first observe that the second component of the array 
-(oo^°(W))"'[ - ooA\W)dkW + »o,(W, V(i)$)] is eq ual to - v'^dkP - i}oo{v,p)dkvK It thus 
follows directly from considering the second component of (|10.3.15p that 

(10.3.17) dtP ^ ~v''dkP-£loc{v,p)9kv'' + 0^-i(c-2; w, V^^^W, $, V^^^*) 

+ 0""i(c-i;W, V(i)W,c-i$,c-iD$) n 0""i(c-2; W, V^^'W, $, 1)$). 

Therefore, since dtp - dtP = {e'^'^l'' - l)dtP ~ 4(c-2aj$)e-4*/c'p, we use Lemma [ES (pXS]) . 
(|6.3.10p . (|6.3.1ip . Lemma (6^ and (|10.3.17p to conclude that 

(10.3.18) dtP^~v''dkp~£loo{v,p)dkv'' + 0^-\c-2;W,v(i)W,$,V«$) 

□ 

Lemma 10-6. There exists a constant C(id) > such that 

(10.3.19) IIUo)llff" ^c-2c(?;d) 

(10.3.20) ||(A - ^2)^0 - 9J(0)||h"-i < c-^C{id). 

Proof. The estimate (|10.3.19p follows from the estimate (|10.3.7p for l{t) at i = O and ()10.3.10p in 
the case m — 0. 

To obtain the estimate (|10.3.20p . first recall that according to the assumption (|8.1.7p and the 
chain rule, we have that 

(10.3.21) (47rG)-i(K2 - A)io = dk{mociv,p)v'') 

= —Q^{f],p)v^dkf]+—^{f],p)v^dkP + 'iy\oo{f],p)dkV^. 

Furthermore, by Lemma [6-21 (110. 3. 8p at < = 0, (110.3. lip in the case m = 0, the chain rule, (|4.1.ip . 
definition (|4.2.6p . (|10.3.18p . and (|3.1.16p in the case c = oo, we have that 

(10.3.22) 

{47TG)-'dtm = ~^{v,p)v''dkV - ^{v,p)v'dkp - 9^00(^,^)9^^^^= + 0''-\c-^■,^d). 
orj op 

The estimate (|10.3.20p now follows from (|10.3.2ip and (|10.3.22p . 

□ 

Lemma 10-7. 

(10.3.23) III d^r, |||h«-2,.. Ill d^P ||k«-2,, < C{id; A,, A^, L,, L'^, L^) L^itd; A,, A2, L,, L'^, L^). 
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Proof. To obtain the bound for d^p, differentiate each side of the expression (|10.3.18p with respect 
to t, and then apply Lemma 16-31 to conclude that 

(10.3.24) dfp^-dt[v''dkP + £lcoiv,P)dkv'']+(Sc, 

where &c G J^-2(c-i; W, DW, Vt^'^tW, V(i)$, c-ift*, c-iV^i^ft^, c'l^^^). We now use 
Lemma lS^ the induction hypotheses, the previously estabhshed bounds (|10.3.14p on 1 1 1 9tW 1 1 
and III dt~W |||hn-i and the definition of I^~^(- • ■ ) to conclude the estimate (|10.3.23p for 

II |||ffiV-2_^ . 

The estimate for d^r] is similar, and in fact much simpler: use equation (|4.1.ip to solve for 9fr/, 
and then differentiate with respect to t and reason as above. □ 

Lemma 10-8. 

(10.3.25) |||/|||^«,, <C(zd;Ai,A2) 

(10.3.26) III dtl |||//«~i,r < C(zd;Ai,A2,Li,4) 

(10.3.27) III dfl \\\H«-2,r < C(zd;Ai,A2,Li,4,L3,i4). 

Proof. To prove ()10.3.25p , we first consider the formula for I given in ()10.3.7I1 + ()10.3.10p . By Lemma 
lOand p0.3.6p . we have that ||| $Koo(?7,p) - $Koo(?7,p) |||ff",T < III ^ooiViP) - ^oo{v,p) \\\ h«,t + 
II ^oo{ti,p) ~ 9\oo{f i,p) |||ff«_^< C(id; Ai, A2). To est imate | || \ \\ H«,r, w here is from (|10.3.7p . 
simply use (I10.3.10p in the case to = together with (|10.3.2[ ) and (|10.3.6p . The proofs of (|10.3.26p 
and ()10.3.27p follow similarly from the expressions (jl0.3.8p . (|10.3.9p . (|lG.3.1ip in the case m = 1, 
and (|10.3.12p in the case m = 1, together with Lemma [6- 2 1 and the bounds supplied by the induction 
hypotheses. Lemma [10-41 and Lemma [10-71 □ 

Lemma 10-9. 

(10.3.28) c"^ III dt<i> IWh^^t < 1/2 + r ■C(i<i;Ai,A2,Li,4) = 4/2 -I- r • C(id; Ai, A2, Li, 4). 

Proof (110.3. 28P follows from definition p0.2.14p . Lemma [1031 inequality (|10.3.26p of Lemma [TOISl 
and inequality (jA.4p of Lemma IA-21 □ 

Lemma 10-10. 

(10.3.29) III 9t$ \\\HN^r<C{id;Ai,A2,Li,L'2)+T-C{id;Ai,A2,Li,L'2,L3,L4) 

= L2{id; Ai, A2, Li, ^2)72 + r • C{id; Ai, A2, Li, 4, L3, L4) 

(10.3.30) III (9^$ \\\HN-i^r < 1/2 + T • C{id;Ai,A2,Li,L'2,L3,Li) 

=^Li/2 + T ■ C{id- Ai, A2, Li, 4, L3, ^4). 

Proof. The estimate (|10.3.29p follows from Lemma [TO^ inequalities (|10.3.26p and (|10.3.27p of 
Lemma 110-81 and inequality (|A.25P of Proposition IA-51 The estimate (|10.3.30p follows from defi- 
nition p0.2.15p , Lemma 110-61 inequality (|10.3.27p of Lemma 110-81 and inequality (|A.6p of Lemma 
E31 □ 

Remark 10.7. Literestingly, (|10.3.29p provides a bound for ||| 9f$ |||_ff".T that depends in part on 
the bound ||| c^^(9f$ |||ff« L'2. 

Lemma 10-11. 

(10.3.31) III $ |||?,«+i,,< + r • C{id- Ai, A2, L2) + r2 • C{id- Ai, A2, L^.L'^^L^, U). 

Proof. Inequality (|10.3.3ip follows from definition (|10.2.13p . (|10.3.25p . (|10.3.29p . and inequality 
(|I30)) of Corollary EH □ 

Lemma 10-12. Let he the energy current (j7.1.ip for the variation W defined by the BGS 
V, and let h {f, g, ■ ■ ■ , ft.*-^-*), where f,g,--- , /i*-'^^ are the inhomogeneous terms from the EOV^ 
satisfied by W that are defined in (110. 2. 4p - ()10.2.6p and that also appear in the expression (|7.3.ip 
for the divergence of'-'^^d- Then on [0,TJ"°^), 

(10.3.32) \\d^{^''>f)\\L^ < C(zd; Ai, A2,ii,4) • [|lW||i. + ||W|U. |lb|U.] . 
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Proof. We separate the terms on the right-hand side of (|7.3.1[) into two types: those that de- 
pend quadratically on the variations, and those that depend hnearly on the variations. We first 
bound (for all large c) the norm of the terms that depend quadratically on the variations by 
C{id; Ai, A2, Li, L'2) ■ ||W||^2- This follows directly from the fact that the coefficients of the qua- 
dratic variation terms can be bounded in L°° by C{id; Ai, A2, Li, Lj). Such an L°° bound may be 
obtained by combining Remark 16.71 Lemma 16-51 in the case m = 1, Remark 16.101 the induction 
hypotheses, (|10.3.14p . and Sobolev embedding. 

We similarly bound the norm of the terms that depend linearly on the variations by 
C{id; Ai, A2) • ||W||i2 ||b||i2, but for these terms, we also make use of the Cauchy-Schwarz inequality 
for integrals. □ 

We also state here the following corollary that will be used in the proof of Theorem ID 

Corollary 10-13. Let V e C\[0,T] x M^) n C^{[0,T],H^ ) n C\[0,Ta], H^'^), and assume that 

V([0,T] X M3) c J?, where K is defined m p0.2.26p . Let W be a solution to the EOV^ (ISlaS]) - 
(|5.0.10p defined by the BGS W with inhomogeneous terms h — {f, g, ■ ■ ■ , ft.^"^-*), where W denotes the 
first 5 components ofV. Let be the energy current (|7.1.2p for the variation W defined by the 

BGSW. Then on [0,T], 

(10.3.33) < C{A;\\\ W \\\l^,t,\\\ dtW \\\l^,t) ■ [\m\h + II^IU^ HblU^] . 

Proof. We do not give any details since CoroUay 110-131 can proved by arguing as we did in our 
proof of Lemma [10-121 In fact, the proof of CoroUav 110-131 is simpler: c does not enter into the 
estimates. □ 

Lemma 10-14. 

(10.3.34) III W ||U«^,< [Ai/2 + T-C(zd;Ai,A2,Li,L^)] • exp(T • C(zd; Ai, A2,Li,4)). 

Proof. Our proof of Lemma 1 1 0- 1 41 follows from a Gronwall estimate in the norm of the variation 
W defined in (|10.2.2p . Rather than directly estimating the i?^ norm of W, we instead estimate 

the L^ norm of '■'^-'tjlj, where '■'^-'tjs is the energy current for the variation (?sW defined by the BGS 
V. This is favorable because of property (I7.2.ip and because by ()7.3.ip . the divergence of ("^^^ is 
lower order in W. We follow the method of proof of local existence from [T8]; the only difficulty is 
checking that our estimates are independent of all large c. An important ingredient in our proof is 
showing that for < \d\ < N and t e [0,T™'''^), we have the bound 

(10.3.35) llbslU^ <C(id;Ai,A2,4)(l + ||W||H«), 

where is defined in p0.2.8p . Let us assume (|10.3.35p for the moment; we will provide a proof at 
the end of the proof of the lemma. 

We now let '■^'^^Ss denote the energy current (|7.1.ip for the variation 95W defined by the BGS V, 
and abbreviating 3s ^'^''^5 to ease the notation, we define £(i) > by 

(10.3.36) £2(^)dcf f ^o,^t,s)d's. 

By (j8.3.ip . and the Cauchy-Schwarz inequality for sums, we have that 

(10.3.37) Ca^^ziml^ < < C^ljml^- 
Here, the value Z — Z{id; A2) is given by p0.3.5p . 

Then by LemmaHEEl (!10.3.35p . (|10.3.37p . with C C{id; Ai, A2, Li, L^), we have 

(10.3.38) 2£^£= V / dJ^.d^s<C- V (||95W||i. + jjasWlU. Hb^lU^) 

CLt Im3 

|5|<Af " \S\<N 

< C ■ (||W||^« + ||W||^«) <C + 
We now apply Gronwall's inequality to (|10.3.38p . concluding that 

(10.3.39) £(t) < [£(0) + Ct] ■ cxp(Ct). 
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Using (I10.3.37P again, it follows from (|10.3.39p that 

(10.3.40) \mt)\\H- < {C^ljWmHg^ + Ct) ■ exp(Ci). 

Recalling that W(0) = Wc-("^Wc and taking into account inequality (j9.0.7p . the estimate (|10.3.34p 
now follows. 

It remains to show (|10.3.35p . Our proof is based on the Sobolev-Moser lemmas stated in Appendix 
IB] and the c— independent estimates of Section [HI With the 5 components of the array b defined by 
(|10.2.4p - (|10.2.6p . we first claim that the term c-A°ds{{cA°y^h) from (|10.2.8p satisfies 

(10.3.41) ||,yi°95((c^")"'b)|U2 < C(*d;Ai,A2,L'2), 

where, as is true throughout this section, all of the estimates we derive are valid on [0, J"™"^). Because 
(16.3. 14p and the induction hypotheses together imply that Hc.A^CVV, $) ||l°o < C{id; Ai, A2), it suffices 
to control the norm of {cA'^)~^h. Then by the induction hypotheses, (|6.3.14p . Proposition IB-21 
and Remark IB. 11 with (c^l^CVV, $))^^ playing the role of F in the proposition and b playing the 
role of G, we have that 

(10.3.42) \\{,A")-^h\\HN <C{id;Ai,A2)\\h\\HN. 
To estimate ||bl|jifw, we first split the array b into two arrays: 

(10.3.43) b = «8c(W, $, £1$) + J^iid, W, $), 

where *8c is defined in Lemma 16-81 and the 5-component array Jc comprises the terms from (jl0.2.4p 
- (|10.2.6p containing at least one factor of the smoothed initial data. By Lemma [6-21 Lemma [6-51 
and Remark l6.101 we have that 

(10.3.44) 3c&T^{id,W,^), 
and from (|10.3.44p and the induction hypotheses, it follows that 

(10.3.45) ||J,(^d,W,$)||^« <C(id;Ai,A2). 

Furthermore, by (|6.3.18p in the case m — 1 and the induction hypotheses, we have that 

(10.3.46) ||<Bc(W,$,D$)||ff« <C(id;Ai,A2,4). 
Combining (|lG.3.43p . (|10.3.45p and (|10.3.46p . we have that 

(10.3.47) ||b||^« <C(zd;Ai,A2,4). 

Now (|10.3.42p and (|10.3.47p together imply (|10.3.4ip . 
We next claim that the kg terms (|10.2.9p satisfy 

(10.3.48) ||ks|U2 <C(^d;Al,A2)||W||^«. 

Since \\cA"{W,^)\\l-^ < C(id; Ai, A2), to prove (|10.3.48p . it suffices to control the norm of 
{cA°)-\A''dk{dsW) - ds (^{cA°)-^cA''dkW^ ■ By the induction hypotheses, (|6.3.14p . Proposition 

and RemarkEl with {cA°)-\A'' = (^{cA")-\A''^ (W, $) playing the role of F in the propo- 
sition, and 9feW playing the role of G, we have (for < |a| < N) that 

(10.3.49) \\{,A'')-\A''ds{dkW) ~ ds{{cAY'cA''dkW)\\L2 < C(zd; Ai, A2)I|V(i)WI|^«-i, 

from which p0.3.48p readily follows. This concludes the proof of (|10.3.35p . and therefore also the 
proof of Lemma 110-141 □ 

11. The Non-relativistic Limit of the EN^ System 

In this section, we state and prove our main theorem regarding the non-relativistic limit of the 
EN^ system. Before stating our main theorem, we first state and prove a corollary of Theorem [2] 
that will be used in the proof of Theorem 21 and we also briefly discuss local existence for the EP^ 
system. 
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11.1. EN|; well-approximates EP^ for large c. The following corollary shows that for large c, 
solutions to the ENJ; system are "almost" solutions to the EP^ system. 

Corollary 11-1. For all large c, the solutions V = (W, $,£>$) to the EN^ system (|4.1.ip - (|4.1.8p 

furnished by Theorem [D satisfy 

(11.1.1) oc^^(W)a,,W = «Boc(W, V(l)$) + (£le 

(11.1.2) A($ - $,) - K^{^ - $,) = A7rG[m^iv,p) - m^ifi,p)] + €2„ 
where 

(11.1.3) III €lc \\\H^-i,T<c-^Ciid-Ai,A2,L2) 

(11.1.4) III £2, \\\HN-i^T<c-'Ciid;Ai,A2,U), 
and T is from Theorem 

Remark 11.1. Note that the left-hand side of (jll.l.ip involves d^W rather than d^W. 

Proof. The estimate ()11.1.3p follows from multiplying each side of (|10.3.16p by oo-A^CW), and then 
combining Proposition IB-2| Remark IB.li (jl0.2.1d|l , and the induction hypotheses from Section 
110.3.11 which are valid on [0, T]. Similarly, the estimate (|11.1.4p follows from the fact that 
A($ - i>c) - - $^) ^ c -^d^^ + l, where / is given by (|10.3.13p . together with p0.3.7p . (|10.3.10p 
in the case m — 0, (jl0.2.1f|) . and the induction hypotheses. □ 

11.2. Local Existence for EP^. 

In this section, we briefly discuss local existence for the EP^ system. 

Theorem 3. (Local Existence for EP^) Let Voc denote the initial data (|8.1.ip for the EP^ 
system (|4.2.ip - (|4.2.6p that are subject to the conditions described in Section\^ Assume further that 
the equation of state is "physical" as described in sections \3.1\ and 1 6. S\ Then there exists aT > 
such that (|4.2.ip - (|4.2.6p has a unique classical solution VooC^jS) on [0, T] x M.^ of the form Voo — 
iVocPoo, ■ ■ ■ $00) TOi/i Voo(0, s) = Voo(s). The solution satisfies Vqc{[0,T] x M^) C .ft, where the 
compact set J? is defined in (|10.2.26p . Furthermore, G C°{[0,T],H^ ) n C^{[0,T], H^'^) and 
<^>eC^{[0,T],H^+^)nCH[Q,T],H^JnC\[Q,T],H^-^). 

Proof. Theorem [3] can be proved by an iteration scheme based on the method of energy currents: 
energy currents (00) 3 can be used to control HWooCi) H^fjy , while || $00 (Oil can be controlled using 

an easy estimate on the operator A — k^. These methods are employed in the proof of Theorem |4] 
below, so we don't provide a proof here. Similar techniques are used by Makino in |12j . We remark 
that these methods apply in particular to the system studied by Kiessling (as described in Section 
1121) in [TU]. □ 

11.3. Statement and Proof of the Main Theorem. 

Theorem 4. (The Non-relativistic Limit of EN);) Let Voc denote initial data (|8.1.ip for 

the EPf^ system (j4.2.ip - (j4.2.6p that are subject to the conditions described in Section O Let 
Vc denote the corresponding initial data (jS.l.Sp for the EN"^ system ()4.1.1I) - (|4.1.8p constructed 
from Voo «s described in Section \^ and assume that the c— indexed equation of state satisfies 
the hypotheses (|6.3.ip and (j6.3.2p and is "physical" as described in sections \3.1\ and \6.3\ Let 
Voo = {'r]oo,Poo,vl^, • ■ • , 93$oo) (Vc = (?7c,Pc, ■ ■ ■ , 83^0)) denote the solution to the EP^ (L^N^) 
system launched by Voo (Vc) as furnished by Theorem\^ (Theorem\^ . By Theorem\^and Theorem 
we may assume that for all large c, Voo and Vc exist on a common spacetime slab [0, T] x R'^, 
where T is the minimum of the two times from the conclusions of the theorems. Let Woo md Wc 
denote the first 5 components of Voo and Vc respectively. Then there exists a constant C > such 
that 

(11.3.1) \\\Woo -Wc\\\h«-kt < -C 

(11-3.2) III ($oo - *oo) - ($c - $c) |||ff« + i,T <C~^-C 

(11.3.3) lim |$oo-$c| =0, 

C— i-CX3 
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where the constants $00 I'rid $c are defined through the initial data by (jS.l.Sp and (|8.1.12p respec- 
tively. 

Remark 11.2. (jll.3.ip . (|11.3.2p . pi.3.3p . and Sobolev embedding imply (recall iV > 4) that 

Woo uniformly and $c — > •I'cx; uniformly on [0, T] x M'^ as c — > 00. Furthermore, the interpolation 

estimate (|B.36I) . together with the uniform bound ||| Wc |||//« C that follows from combining 

(I6.3.9p . pp. 2. lap , and (|10.2.1bp . collectively imply that limc^oo ||| Woo — Wc \\\hn' t= for any 
N' < N. The reason that we cannot use our argument to obtain the norm on the left-hand 
side of (|11.3.ip instead of the H^~^ norm is that the expression (|11.3.8p for b already involves one 
derivative of W, and therefore can only be controlled in the H^^^ norm. 

Proof. Throughout the proof, we refer to the constants Ai, A2, etc., from the conclusion of Theorem 
[21 but we typically suppress the dependence of the running constants C(- • • ) on N, k,02,02, K, 
and To further ease the notation, we drop the subscripts c from the solution Vc and its first 5 
components W^,, setting V Vc, W W^,, etc. We then define with the aid of (|8.1.14p 

(11.3.4) W''=^'Woo-W 

(11.3.5) $ = ($00 - $00) - ($ - *c) - (*oo - *oo) - (* - $c). 

Our proof of Theorem [4] is similar to our proof of Lemma I10-14| we use energy currents and ele- 
mentary harmonic analysis (i.e. Lemma rA-4p to obtain a Gronwall estimate for the H^^^ norm of 
the variation W defined in (|11.3.4p . It will also follow from our proof that the H^^^ norm of $ is 
controlled in terms of ||W||jLfN-i plus a small remainder. We remark that all of the estimates in this 
proof are valid on the interval [0,T], where T is as in the statement of the theorem. 

From definitions (|11.3.4p and (|11.3.5p . it follows that IV, $ are solutions to the following EOV^ 
defined by the BGS Woo : 

(11.3.6) oo.A^(Woo)a^W = b 

(11.3.7) (A-k2)$==;, 
where 

(11.3.8) b = 'Boo(Woo, VW$oo) - »co(W, V«$) + [oo.A^(W) - oo^^(Woo)]a^W - €lc, 

(11.3.9) Z = 47rG[5noo(ryoo,Poo)-$noo(??,p)] - €2^, 

*8oo is defined in Lemma [6-81 and £lc, 2;2c are defined in CoroUarv lll-ll Note that the definition 
of I in (|11.3.9p differs from the definition of I (|10.3.13p used in the proof of Corollary 111-11 By 
comparing (|8.1.ip and (|8.1.8p . we see that the initial condition satisfied by W is 

(11.3.10) W(t = 0) = O. 
Differentiating equation (|11.3.6p with 9^, have that 

(11.3.11) oo^^(Woo)9^(9alV) =b5, 

where (suppressing the dependence of oo-A'^{-) on Woo for v — 2,3) 

(11.3.12) bs = ooA^ds ((oo.A°)-ib) + ka 
and 

(11.3.13) ks = oo^°[(oo.A°)-ioo.A'=95(9feW) --as((oo.A°)-ioo^'=afcW)]. 

As an intermediate step, we will show that for < jo?! < iV — 1 and t E [0, T], we have that 

(11.3.14) llbslU- < C{id;\\\ Woo |||^^«^,T,Ai,A2,ii,L2,L4) • {\\W\\h"-^ + cr') . 

Let us assume ()11.3.14p for the moment and proceed as in Lemma 110-141 we let '^°°^ds denote 
the energy current (j7.1.2p for dsW defined by the BGS Woo, and define £(<) > by 

(11.3.15) £2(0= E / 3%it,s)d\ 

\a\<N-l''^^ 
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where we have dropped the superscript (cxd) on 2 to ease the notation. By (|8.3.ip and the Cauchy- 
Schwarz inequahty for sums, we have that 

(11-3.16) Co,,zm?HN-^ < < C^lJW\\'H^-r- 

Here, Z =\\\ ^^o |||l~,t ■ Then by Corollary [TOUS] + Sobolev embedding, (|11.3.14p . and (|11.3.16p . 
with C = C(id; \\\ Woo |||_f/« .tj III dtWoo |||ff"-i,Ti ^ii ^2, ^i, -^2, L4j, we have that 

(11-3.17) 2£|£= I «c?3s<C- (ll^s^ll^. + liasWIU.IIbslU.) 

|S|<Ar-l \a\<N-l 

< C - ||W||^„_i + c-^C - ||W||^«-i < C - £2 + c-^C - £. 

Taking into account (|11.3.10p . which implies that £(0) = 0, we apply Gronwall's inequality to 
(|11.3.17p . concluding that for te [0, T], 

(11.3.18) £{t) <c-^C -t-expiC -t). 

From (|11.3.16p and pi.3.18p . it follows that 

(11-3.19) III IV \\\h'^-\t^ c~^C ■ T ■ exp(r - C), 

which implies (|11.3.ip . 

We now return to the proof of (|11.3.14p . To prove (|11.3.14p . we show only that the following 
bound holds for t G [0,T], where for the remainder of this section, 
C = C[id] III Woo |||//«-i T'^i' ^2, -^1,-^2,^4) : 

(11.3.20) ||b||ffiv-i < C ■ ||W||ffiv-i + c~V. 

The remaining details, which we leave up to the reader, then follow as in the proof of Lemma [10-141 
By (|10.3.6p . which is vahd for t = T, and by (|R32|, we have that 

(11.3.21) ||5noo(?7oo,Poo) - 5noo(?7,p)||H"-i < C ■ ||W|l^«-i, 

and combining pi.l.4p . (|11.3.7p . (|11.3.9p . (|11.3.2ip . and Lemma [Ol it follows that 

(11.3.22) ||^||h"+i < C ■ \\1\\h'^-i < C ■ ||W||^iv^i + c-^C. 
Similarly, taking into account (|11.3.22p . we have that 

(11.3.23) ||Q5oo(Woo, V(i)$oo) - »oo(W, v(i)$)||h«-i < C ■ (||W||^«-i + I|VW*I|h«-i) 

<C- ||W||h"-i +c-^c. 

Finally, by (|10.3.6p and (|10.3.14p . which are both valid for r = T, by ((R30)) . and by ((R32)) . we 
have that 

(11.3.24) ||[oo.A'^(W)-oo^^(Woo)]a^W||^„_, <C.||W||^«-i. 

Inequahty (|11.3.20p now follows from (|11.1.3p . pi.3.8p . (|11.3.23p . and (|11.3.24p . The estimate 
(I11.3.2P then follows from (|11.3.5p . (|11.3.19p . and (|11.3.22p . while (|11.3.3p is merely a restatement 
of (|8.1.13P . □ 
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A. Inhomogeneous Linear Klein-Gordon Estimates 



In this Appendix, we collect together some standard energy estimates for the linear Klein-Gordon 
equation with an inhomogeneous term. We provide some proofs for convenience. 



Lemma A-1. Let I G C°([0, T], and «'o(s) £ H^(R^), where N e N. Then there is a 

unique solution $(t,s) : R x — R to the equation 

(A.l) -c-^df^ + A$ - = Z 

with initial data $(0,s) = 0, 9f$(0, s) = 5'o(s). The solution has the property 

$ G c°([o,T],i7^+i(R3)) ncH[o,r],i/^(R3)). 

Proof. This is a standard result; consult 17 \ for a proof. □ 
Lemma A-2. Assume the hypotheses of Lemma \A-l\ Assume further that 

dtl e C°([0,T],i7^-i(R3)) anddfl G C°{[0,T], H'^-^{R^)). Then there exists a constant Co{k) > 
such that 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 
(A.8) 



'i' |||/i""+i,T< Co(k) • (c "^II'I'oIIhn -f cT III Z lll/i-N t ) 

dt^ ||U«,T< Coin) ■ {W-^oWh^+c^T III I IWhn^t ) 

|||ff«,T< Coin) ■ (ll*o||//« +c||/(0)||ff«-i -KcT ||| dtl \\\h^-i,t ) 
\\\hn-kt< Coin) ■ {c\\i'o\\HN +c'\\liO)\\HN-i+c'T \\\ dtl ||k«-i,T ) 
|||h«-i,t 

< Coin) ■ {c'WmWH^-i + c\\iA - n^)i'o - dtmWn^-^ + cT \\\ dfl \\\h«-2^t) 
Of^ 1 1 \h^~'^ t 



<CoiK)-ic'\\liO)\\H^ 



'-WiA - _ a,/(0)||^«-. + c^T III d^l \\\h 



N-2 7^ 



Proof. Because V'^'^'<i' is a solution to the Klein-Gordon equation 

—c^'^dtiy^'^^^) + A(V^''')$) — K^(V(''')$) = V^'^H, we will use standard energy estimates for the 
linear Klein-Gordon equation to estimate ||| $ |||jLfN+i rp . Thus, for < fc < A^, we define E^it) > 
by 

dcf 



(A.9) 



Elit) = ||«:vw$(t)||i2 + iivc^'+i^^iii. + ||c-ivwat$(i)|| 



2 
L2- 



We now multiply each side the equation satisfied by V'-'"'^$ by —V^'^^dt^, integrate by parts over 
R'^, and use Holder's inequality to arrive at the following chain of inequalities: 



(A.IO) 



E,it)^^E,it) = (Elit)) = ( - vWa.<i>) • (V^'^h) ds 

< ||vWat4>(t)|u.||vW/(t)|U2, 



where ( - V'-'^'^dt^) ■ (V'-'^H) denotes the array-valued quantity formed by taking the component by 
component product of the two arrays — V'''^'9t$ and V'-'"'/. 
If we now define Eit) > by 



AT 



(A.ll) E'it)'='iJ2EUt) 



it follows from (jA.lOp and the Cauchy-Schwarz inequality for sums that 
(A.12) 
and so 
(A.13) 

Integrating (jA.13p over time, we have the following inequality, valid for t G [0, T] 
(A.14) Eit)<EiO) + ct\\\l\\\HN^T- 



mjEit) = —iE'it)) < ||a,<i>||H«||/(t)||^« < cEit)\\lit)U., 



-Eit)<c\\lit)\\H^. 
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From the definition of E{t) and the initial condition $ = 0, we have that 

(A.15) \\^{t)\\HN+^ <C{K)E{t) 

(A.16) WdtmU- < cE{t) 

(A.17) £;(0) =c-i*o||H«- 

Combining (|A.14p . (|A.15p . (|A.16p . and (|A.17p . and taking the sup over t e [0,T] proves (fO)) 
and (|A.3p . 

To prove (|A.4p - (|A.7p . we differentiate the Klein-Gordon equation with respect to t (twice to 
prove (|A.6P and (|A.7P ) and argue as above,taking into account the initial conditions 

(A.18) 5(2$(0) = -c^l{0) 

(A.19) d!^0) = [(A - K^)i'o - dtm] . 

□ 

Corollary A-3. Assume the hypotheses of Lemma lA-Si and let Co(k) be the constant appearing in 
the conclusions of the lemma. Then 

(A.20) III $ \\\jiN+i^T< {Co{k.))^ ■ {c-^\\i'a\\lN + 2 • T- ||| 9*$ \\\hn,t ■ \\\ I \\\h«,t )• 

Proof. Inequality (|A.12p gi ves that ^jt{E'^{t)) < \\dt^H"\\li't)\\H'^ ■ Taking into account (|A.15P and 
(|A.17p . the proof of (IA.20p easily follows. □ 

Lemma A-4. Let TV e N, and 3 e H^-^{R^). Suppose that $ e H^+\R^) and that A$-k2$ = 3. 

(A.21) ||$||ff« + i(R3) < C{N, At)||a||jy«-l(R3). 

Proof. For use in this argument, we define the Fourier transform through its action on integrable 
functions F by F{£^) J^^ F(s)e~^'^*^'® ds. The following chain of inequalities uses standard results 
from Fourier analysis, including Plancherel's theorem: 

(A.22) \mi^ < c\\{i + < c{k) [ («2 + |2^en'i%)prf'e 

and this proves (|A.2ip in the case = 1. To estimate norms of the fc*'* order derivatives of $ for 
fc > 1, we differentiate the equation k times to arrive at the equation A(V'^ — k^(V^'^'<I>) = V^'^-'U, 
and argue as above to conclude that 

(A.23) l|V('=^$||^2 <C{K)\\^^''b\\l,. 

Now we add the estimate (|A.22I) to the estimates (|A.23I) for 1 < fc < - 1 to conclude (|A.2ip . □ 

Remark A.l. The hypothesis $ G iJ^+^(R'^) does not follow from the remaining assumptions. 
For example, consider g{x) = e"". Then g - -^g e i^(M), but g ^ i?^(R). 

Proposition A-5. Assume the hypotheses of Lemma \A-l\ Assume further that 

I e C0([0,T],i7^(M3)), dtl e C°{[0,T],H^-\R^)), and d^l £ C°{[0,T], H^-^{R^)). Then 

(A.24) III $ \\\h'-'+\t < C{N,k) 

■ (c-i*o||ff« + ||/(0)||h«-i+ III I \\\h«~\t + T III dtl \\\h«~\t ) 

and 

(A.25) III 9t$ lllfl-iv^T < ^(A^,^) 

• (c||/(0)||ff«-i + ||(A - _ at/(0)||^«-2+ III dtl \\\h«-^^t +T III d^l \\\h«-2^t )• 

Proof. Define 3 I + c^^df^ and observe that $ is a solution to 
(A.26) A$ - a. 
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By inequality (jA.Sp of Lemma lA-21 Lemma rA-41 and the triangle inequality, we have that 
(A.27) III $ |||ffiv+i,T< C{N,k) III l + c-^d^^ \\\h'^~\t 

< C{N, K){c-^\\i>o\\H^ + ||/(0)||ff«-i+ III I \\\hm-i^t + T III dtl |||ff«-i,T ), 

which proves (|A.24[) . 

Because 9t<i> satisfies the equation —c^^d^^ + A(9t<l') — K^(9t<i>) — dtl, we may use a similar 
argument to prove (|A.25|) : we leave the simple modification, which makes use of (|A.7|) . up to the 
reader. □ 

B. SOBOLEV-MOSER ESTIMATES 

In this Appendix, we use notation that is as consistent as possible with our use of notation in the 
body of the paper. To conserve space, we refer the reader to the literature instead of providing proofs: 
propositions IB-21 and IB-41 are similar to propositions proved in chapter 6 of [9] , while Proposition 
IB-51 is proved in [TT]. The corollaries and remarks below are straightforward extensions of the 
propositions. With the exception of Proposition IB-61 which is a standard Sobolev interpolation 
inequality, the proofs of the propositions given in the literature are commonly based on the following 
version of the Gagliardo-Nirenberg inequality [13], together with repeated use of Holder's inequality 
and/or Sobolev embedding: 

Lemma B-1. Ifi,kGN with < i < fc, and V is a scalar-valued or array-valued function on M.'^ 
satisfying V G L°°(K'') and || V^''^ V||^2(Rd) < oo, then 

(B.28) ||V«V|L../. < C(fc)||v||^-*||VWV|||,. 



Proposition B-2. Let K C W be a compact set, and let j, d G N with j > ^. Let V : R'^ ^ W 

be an element of (U.'^), and assume that V(R'') C K. Let F G Cl{K) be a q x q matrix-valued 
function, and let G G H^iW^) be a q x q (q x 1) matrix-valued (array-valued) function. Then the 
q X q {q X 1) matrix-valued (array-valued) function (Fo V)G is an element o/iJ-' (R'*) and 

(B.29) \\{F o V)G||^.(R.) < CU, d)\F\,,K{l + ||V|P^,(«,))||G||^,(r.). 

Corollary B-3. Assume the hypotheses of Proposition \B-B with the following changes: V, G G 
C^{[0,T], H^{R'^)). Then the q x q [q x 1) matrix-valued (array-valued) function {F o V)G is an 
element of C°{[0,T], {R"^)). 

Remark B.l. We often make use of a slight modification of Proposition lB-2l in which the assumption 
V G H^(U.'^) is replaced with the assumption V G H^{M.'^), where V G M" is a constant array. Under 
this modified assumption, the conclusion of Proposition [B-21 is modified as follows: 

(B.30) ||(FoV)G||ff. <G(j,d)|FU(l + ||V|p^^)||G||H.. 

A similar modification can be made to Corollarv lB-31 

Proposition B-4. Let K C M" be a compact, convex set, and let j,d eN with J > |. Let F £ Cl [K) 
be a scalar or array-valued function. Let V, V : R*^ ^ R", and assume that V, V G (W^). Assume 
further that V(R''), V{R'^) C K. Then FoV-FoVg W {R'^) and 

(B.31) ||F o V - F oV\\h, < G(j, d, \\V\\h, , ||V||^,)|^^|,+i,k||V - V\\h, ■ 

Remark B.2. As in Remark lB.li we may replace the hypotheses V, V G H^(R.'^) from Proposition 
IB-41 with the hypotheses V, V G H^{M.'^), in which case the conclusion of the proposition is: 

(B.32) ||(FoV)-(FoV)||h. <GO;d,||V||^,,||V||^,)|^^|,+i,K||V-V||^.. 

V V 

Furthermore, a careful analysis of the special case V = V, where V G iiT is a constant array, gives 
the bound 

(B.33) ||FoV-FoV||h. <G(j- d)|aF/aV|,_i,K(l + ||V|p^-^i)(||V||^^), 

in which we require less regularity of F than we do in the general case. 
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Proposition B-5. Assume the hypotheses of Provosition [B^ with the following two changes: 

(1) Assume j > | + 1. 

(2) Assume that G £ W-^iR"^). 

Let fc e N with 1 < k < j, and let a be a spatial derivative multi-index with \d\ — k. Then 

\\d3{{FoV)G)-{FoV)dsGh2 
(B.34) < Cad)|9F/aV|,_i,K(||V||ff. + \\vr^,)\\G\\H.-^. 

Remark B.3. As in Remark lB.ll we may replace the assumption V £ H^{M.'^) in Proposition IB-51 
with the assumption V G Hl^{M.'^), where V is a constant array, in which case we obtain 

\\ds{{F oV)G) - {F oV)dsG\\L2 
(B.35) < C{j,d)\dF/dV\,^,.K{\\V\\^, + \\V\^)\\G\\h.-^. 

Proposition B-6. Let N',N eR be such that < N' < N, and assume that F e i/^(M'0- Then 
(B.36) < C(7V',d)||F||i-^'/^||F||2:/^. 
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